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Time reversal symmetric topological superconductors in three spatial dimensions carry gapless 
surface Majorana fermions. They are robust against any time reversal symmetric single-body per¬ 
turbation weaker than the bulk energy gap. We mimic the massless surface Majorana’s by coupled 
wire models in two spatial dimensions. We introduce explicit many-body interwire interactions that 
preserve time reversal symmetry and give energy gaps to all low energy degrees of freedom. We show 
the gapped models generically carry non-trivial topological order and support anyonic excitations. 


I. INTRODUCTION 

Topological superconductors (TSC) are electronic 
phases of matter with finite excitation energy gaps that 
are not continuously connected to a conventional BCS 
5-wave superconductor. In particular BCS superconduc¬ 
tors in three dimensions can have non-trivial topologies 
protected by time reversal symmetry.^^^^^ There is a bulk 
integral quantity N of the mean-field system, known as 
chirality, that cannot change upon any adiabatic evo¬ 
lution unless the energy gap is closed or time reversal 
symmetry (TRS) is broken. TSC also exhibits unique 
physical signature along its surface. Despite there is a 
bulk energy gap, the surface of a TSC hosts N gapless 
Majorana (real) fermion modes that are robust, in the 
single-body mean-field framework, to all symmetry and 
bulk gap preserving perturbations. The superfiuid ^He- 
d^lEland perhaps superconducting Cua;Bi 2 Se 3 EEIare can¬ 
didates of TSC. 

The Z classification of TSC - or class Dill band the¬ 
ories according to the Altland-Zirnbauer classificatiorP^ 
- relies heavily on the single-body BCS description of 
the electronic structure. It has recently been shown that 
under strong many-body interaction, the surface state of 
sixteen copies of a TSC can be gapped without breaking 
time reversal symmetry or introducing surface topologi¬ 
cal order. This reduces the integer classification of TSC 
into Zi6.^^^®^This suggests the many-body extension al¬ 
lows a continuous path that connects sixteen copies of 
a TSC to a trivial 5-wave superconductor in three di¬ 
mensions without breaking symmetry or closing the bulk 
gap. In fact, the surface Majorana modes of any TSC can 
be gapped without breaking symmetries. However, there 
would generically be a residue topological order, unless 
A/" is a multiple of 16, that all ows n on-trivial anyonic ex¬ 
citations to live on the surface.^mi^ As a result, these 3D 
bulk systems are still topologically distinct from a trivial 
state. 

Similar phenomena were also seen in topological 
insulators^^J^ in three dimensions and topological 
superconductor^^ in one dimension. Many-body in¬ 
teractions allow the surface Dirac mode of a topologi¬ 
cal insulator to acquire an energy gap without breaking 
time reversal or charge conservation symmetries. How¬ 
ever a non-trivial surface topological order would be left 


behind.l^^M^ This indicates the bulk insulator still car¬ 
ries a non-trivial Z 2 symmetry protected topology (SPT) 
even in the many-body framework. On the other hand, 
the Z classification of time reversal symmetric BDI su¬ 
perconductors in one dimension breaks down to Zg in the 
presence of strong interaction.^^^^^ 

The topological order of a gapped symmetric surface 
of a topological insulator or superconductor was deduced 
mainly using vortex condensation or other topological 
field theory techniques. They do not specify the micro¬ 
scopic many-body surface gapping interactions that give 
rise to these exotic surface states. A pioneer work that 
addressed this issue was done by Fidkowski and Kitaev 
in RefEU where they constructed explicit time reversal 
symmetric 4-fermion interactions that give an energy gap 
to eight boundary Majorana zero modes of a ID TSC. 
Another insightful work was published by Mross, Essin 
and Alicea in RefEU where they mimicked the surface 
Dirac mode of a topological insulator using a coupled 
wire model and wrote down explicit symmetric gapping 
interactions that lead to different gapped or gapless sur¬ 
face states. 

Sliding Luttinger liquid^^^Jt^ and coupled wire 
construction^^ are immensely powerful in building two 
dimensional topological phases. They model 2D systems 
by arrays of coupled ID chains, where interaction effects 
are more controlled and better understood. This theoret¬ 
ical technique has been frequently used in the stu dy of 
fractional Quant um Hall state^^^®Il, anyon modelj^^^, 
spin liquidl^^^, (fr action al) topological insulatorJ^^H^ 
and superconductorJ^^^. 

In this article, we imitate the surface Majorana modes 
of a 3D topological superconductor using a coupled Majo¬ 
rana wire model, construct explicit 4-fermion interactions 
that lead to a finite excitation energy gap, and study the 
residue surface topological order. 


A. Summary of results 

We consider a 2D array of chiral Majorana wires, each 
of which carries N Majorana fermion channels that prop¬ 
agate in a single direction. The chiralities of wires al¬ 
ternate so that adjacent wires counter-propagate and 
Majorana’s can backscatter to their neighbors through 
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electron tunneling (see figure [^. When the interwire 
backscattering is uniform, the 2D system is gapless. In 
the long wavelength continuum limit, the energy spec¬ 
trum is linear in both kx and ky directions and the model 
gives N Major ana cones. 
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FIG. 1. (Left) Coupled wire model ( |2.4| of N ga pless surface 
Majorana cones. ( Righ t) Fractionalization (3.4) and couple 
wires construction ( |3.8[ ) of gapped anomalous and topological 
surface state. 


If the N massless Majorana species decouple, each 
of them is protected by a non-local “antiferromagnetic” 
time reversal symmetry (TRS) that translates all Ma- 
jorana’s to the next wire while reversing the propagat¬ 
ing direction. This symmetry requires uniform interwire 
backscattering and forbids a fermion mass. However, this 
TRS is qualitatively different from the conventional one, 
which is local. For instance a gap can be introduced 
when N is even by inter-species single-body backscatter- 
ings which preserves symmetry.l^ Thus the surface classi¬ 
fication becomes Z2 instead of Z even in the single-body 
framework. Despite the discrepancy, the coupled wire 
model does bare resemblance to the original problem of 
topological superconductor surfaces, especially when the 
number of chiral species N is odd. 

A major result of this work is to construct many-body 
gapping potentials that freeze out all low energy degrees 
of freedom while preserving the non-local time reversal 
symmetry. This is achieved by fractionalizing or hiparti- 
tioning the Majorana channels on each wire into a pair of 
independent sectors without interfering each other. They 
can then be backscattered to adjacent wires in opposite 
directions. When there are even Majorana species, the 
decomposition is obvious as = 2r = r + r and one 
can simply separate the first r Majorana’s from the re¬ 
maining r. The fractionalization in the odd case is more 
involved but are well-known in the conformal field theory 
(CFT) community. Firstly the emergent rotation symme¬ 
try among the fermion species corresponds 

to a chiral so{N) current algebra at level 1, also known 
as an afhne Kac-Moody algebra or Wess-Zumino-Witten 
(WZW) theory, along each wire.^^ Low energy excita¬ 
tions along the chiral wires, referred as primary fields, are 
irreducible representations of the so{N)i algebra. Frac- 
tionalization of the WZW CFT is also known as level- 


rank dualit 3 l^^H^ or conformal embeddin ^^^ l ^^H ^ 

so{N)i ( 1 . 1 ) 


where the two are mutually commuting subalge¬ 

bras of so{N). For example so(9)i can be decomposed 
into the tensor product 50 ( 3)3 G) 50 ( 3)3 as 9 = 3 x 3. 
This splits each wire into a pair of fractional channels 
(see figure [^. For instance the chiral central charge, 
c_ = N/2, which loosely speaking counts the degrees of 
freedom and characterizes the heat currenfP^H^ running 
along each wire, also decomposes so that each Qn channel 
carries c_ = N/A. 

The many-body gapping interaction are given by in¬ 
terwire current-current backscattering (see also figure 


T^int — ^ 




( 1 . 2 ) 


where are the currents operators along wire y. 

All current operators are certain combinations of fermion 
bilinears, and the backscattering interaction therefore 
consists of 4-fermion terms. This Hamiltonian is exactly 
solvable. It preserves the “antiferromagnetic” time re¬ 
versal symmetry and opens up an excitation energy gap. 

The symmetric gapped surface generically carries a 
non-trivial Gn topological order. 


r 5'0(r)i, for N = 2r 

\ S0{3)3 ^6 SO{r)i, for AT = 9 + 2r 


(1.3) 


where both N and r can be extended to negative in¬ 
tegers. It ca n be inferred, using the bulk-boundary 
correspondenc^^^SMUl^ from the (1 + 1)D Qn CFT liv¬ 
ing along the interface that separates the TR symmet¬ 
ric topological gapped domain and a TR br eaking trivial 
gapped domain. The anyon structur^^H^^H^, which en¬ 
codes the quasiparticle types together with their statis¬ 
tics and fusion properties, follows a 32-fold periodicity in 
the sense that Gn — Gn-\- 32 - Moreover, these thirty-two 
topological states exhibit a natural Z32 relative tensor 
product structure, Gn^ ^hGN 2 — Gn-i+N 2 i where cer¬ 
tain set b of non-trivial bosons are condensecP^ under 
the tensor product. 

It is important to clarify at this point that the coupled 
wire construction is a (2 + 1)D model where the inter¬ 
action (1.2) is built out of local bosonic current opera¬ 
tors J. Under this interpretation, the coupled SO{N)i 
wire model is bosonic and Majorana fermions are treated 
as anyonic excitations that carries a quasiparticle string. 
There is a Z 2 gauge degree of freedom that couples to the 
fermions, and there are deconfined Tr-fluxes 

(or he/2e-fluxes), which are anyonic excitations that are 
non-local with fermions. 

When A/" is a multiple of four, the Gn topological or¬ 
der is Abelian with four distinct anyon types 1, 7 /, 5 +, 5_, 
where s± are 7r-fluxes with opposite fermion parities. 
When N is 2 mod 4, the Gn state resembles an Ising 
topological order with anyons 1,7/, a. When N is odd. 
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the topological state has 7 anyon types, 1, <a±, 7±,/5,/, 
and has a structure similar to SO{3)s (or equivalently 
SU (2)6). All these anyon theories contain Tr-fluxes, which 
should be absent on the surface of a fermionic topolog¬ 
ical superconductor. In ReflTTl the surface topological 
order of a A" = 1 fermionic TSC only contains 4 quasi¬ 
particles 1 , 7 ±,/ instead of 7. The additional 7r-fluxes 
in our coupled wire model could become confined by 
re-introducing single-body interwire fermion backscatter- 
ing. In this case, the thirty two bosonic topological states 
reduces down to two fermionic ones, (1) a trivial state 
containing similar to copies of Px + ipy supercon¬ 
ductors when N is even, or (2) a non-trivial fermionic 
SO{3)s state with anyons 1 , 7 ±,/ when N is odd. This 
Z 2 classification, instead of Z 16 , is a natural consequence 
of the “antiferromagnetic” time reversal symmetry. In 
particular there is no reason to expect the result would 
match that of Ref ITT] and [12] when N is even. 

We will introduce the single-body coupled Majorana 
wire model at the beginning of section |TTj A review on 
the so{N)i WZW CFT will be given in section IIA and 
TTbI as well as in appendix Em and 0 In section |inl 


we will construct time reversal symmetry 4-fermion in¬ 
teractions that will open up an excitation energy gap. 
The discussion will be decomposed into the even and 
odd N cases in section [III A| and |IIIB| respectively. In 
the even case, the gapp ing H amiltonian will match the 
0(r) Gross-N eveu mo deP^^^and we will show an energy 
gap in section III A 1 by (partially) bosonizing the prob¬ 
lem. The gapping potential for the odd case will rely on 
a conformal embedding and relat e to t he Zamolodchikov 
and Fateev Ze parafermion CFTT^^^. This will be dis¬ 
cussed and reviewed in section IIIIB 11 IIIIB 21 as well as 


in appendix [P] The symmetric gapping interactions will 
correspond to non-trivial surface topological orders. This 
will be discussed in section [IV| where we will present the 
class of 32-fold periodic topological Gn states. In sec¬ 
tion [V| we will describe alternative gapping interactions 
that would lead to even more possibilities. Lastly, we 
will conclude the article in section ED where we will also 
discuss some possible future exploration. 


II. COUPLE WIRE CONSTRUCTION OF 
SURFACE MAJORANA CONES 


A time reversal symmetric BCS superconductor is de¬ 
scribed by a Bogoliubov - de Gennes (BdG) Hamilto¬ 
nians HBdcO^)' Symmetries require THBdG0^)T~^ = 
-ffBdG(-k) and CiJBdG(k)C“^ = -iJBdG(-k) where T 
and C are the antiunitary time reversal and particle-hole 
operators. When the symmetries square to = 

1, the BdG theory belongs to the symmetry class Dill 
according to the Alt land-Zirnbauer classificatiorP^ and 
theories in three dimensions with finite excitation energy 
gaps are topologically classified by integer^®^. Supercon¬ 
ducting ^He in the B-phas^®^ and certain doped topo¬ 
logical insulatorJ^ were suggested to carry non-trivial 


topologies. 

Topological superconductors host protected gapless 
surface Majorana modes. The simplest version is a single 
Majorana cone, which is the spectrum of a massless two- 
component real fermion 'H± = ^^ where = 

dyTx ± dxTz and the Pauli matrices Tx^Ty^Tz act on the 
surface real fermion -0 = Majorana fermions 

are hermitian and obey the anti-commutation 

relation {'^^(r),-0^/(r')} = 25jj'5{Y — r'). Time rever¬ 
sal switches the components T(<ai'0L + gl 2 '^r)T~^ = 
0 ^ 2 "^L — so that = —1. The sign in the Hamilto¬ 
nian 1-L± determines its chirality. A general surface state 
could consist of multiple copies of Majorana cones with 
different chiralities 

Nr Nl 

n^ = Y^ iValkl^+i’a + Z (2-1) 

a=l 6=1 

Fermions -0^ and with opposite chiralities can anni¬ 
hilate each other by the time reversal symmetric mass 
term imip^Tz'il^i,. Quadratic terms among fermions of 
the same chirality would however either break time re¬ 
versal or only move the gapless Majorana cones away 
from zero momentum without destroying them. The net 
surface chirality N = Nr — Nr is thus a robust topo¬ 
logical signature that distinguishes and characterizes 3D 
bulk topological superconductors. It cannot be altered 
by any time reversal symmetric two-body perturbations 
that are not strong enough to close the bulk excitation 
energy gap. 

Recent theoretical studies suggest many-body interac¬ 
tions can remove these gapless surface degrees of free¬ 
dom. To construct explicit gapping terms, we turn to 
an anisotropic description of surface Majorana fermions 
using an array of coupled fermion wires (see figure . 
The horizontal wires are labeled according to their ver¬ 
tical position y = ... ^ —2, — 1, 0,1, 2 ,... and each carries 
N chiral (real) Majorana fermions 'ij^y = (-0^,...,-0^) 
which propagate only to the right (or left) if y is even 
(resp. odd). The number of flavors N here is going to be 
identified with the net chirality of the surface Majorana 
cone. Time reversal symmetry is non-local in this model 
as it relates fermions on adjacent wires that propagate in 
opposite directions, 

/ N \ N 

^ '7'”' = (-1)"E<7+1- (2-2) 

\a=l / a=l 

Similar to the symmetry of an antiferrormagnet, here 
time reversal on the single-fermion Hilbert space squares 
to a primitive translation up to a sign, = —ty for 
iy the vertical lattice translation y ^ y ^2 that relates 
nearest co-propagating wires. In the many-body Hilbert 
space, 

R = i-lfiy (2.3) 

where (—1)^ is the fermion parity operator whose sign 
depends on the eveness or oddness of fermion number. 
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We mimic N copies of surface Major ana cones by the 
coupled wire Hamiltonian 

oo 

n0= (2.4) 


where the Wcomponent Majorana fermion -0 disperses 
linearly (for small ky) with velocities Ux, Vy along the hor - 
izontal and vertical axes (see figure]^. By applying ( |2.2[ ), 
we see T'HoT~^ = Ho and the coupled wire model is 
therefore time reversal symmetric. Moreover, Ho has 
continuous translation symmetry along x and discrete 
translation al ong y ^ y + 2. The alternating sign in the 
first term of (2.4) specifies the propagating directions of 


the wires. Projecting to the k x = 0 zero modes along 
the wires, the second term in (2.4) effectively becomes 
a ID Kitaev Majorana chairP^ which has a linear spec¬ 
trum for small ky. More explicitly, by using the Nambu 
basis 

the Fourier transform of the Dirac fermion Cy{x) = 
('02y-i(^) + the coupled wire Hamiltonian 

(2.4) can be expressed as Ho = ^k^BdG(^)^k 5 where 

the BdG Hamiltonian is given by 

-f^BdG(k) = 2Vxfca;Ta; + Vy [-Sm kyTy + (1 - COS ky)T^] 

(2.5) 

for —oc < kx < oo and —tt < ky < it. It has a linear 
spectrum near zero energy and momentum as shown in 
figure 



FIG. 2. Th e energy spectrum of the coupled Majorana wire 
model \2.A\ 


We no tice in passing that if the time reversal oper¬ 
ation in (2.2) was defined without the alternating sign 
(—1)^, it would sqaure to a different sign = ^ty i n the 
single-fermion Hilbert space and the vertical term in (2.4) 
would need to be modified into 

order to preserve the symmetry. This would correspond 
to an alternating Majorana chain in the ^-direction, 
where the gapless Majorana cone would be positioned 
dX ky = t: instead of 0 and would still be protected 
by Kramers theorem as = —1. This sce¬ 

nario is actually equivalent and related to the original by 
a gauge transformation ('04y,'04^+2,V^^+i+s) ^ 
(' 04 y,' 04 y+i,-' 04 y+ 2 ,-' 04 y+i+ 3 ), and therefore the sign 


of 7^ is unimportant in this problem. Nevertheless we 
will stick with previous convention defined in (2.2) in the 
following discussions. 

Th e chirality N of the coupled Majorana wire model 
(2.4) is set by the chiral central charge c_ = N/2 
along each wire. This quantity is defined by the dif¬ 
ference of central charges^ between right and left mov¬ 
ing modes, and determines the energy (thermal) cur¬ 
rent It ~ c_ ^ flowing along the wire in low 


temper at urd^^^^. In general, a Majorana wire carry¬ 
ing Nr right moving fermions and Nr left moving ones 
has the kinetic Hamiltonian H = where = 

[JIatp 0 (—ffATr^j ^a; a cts on the (A/'i^+ A/'L)-component real 
fermion -0. In ( |2.4[ ) we consider the simplest case when 
(A/'i^, Nl) = (A^,0) for y even or (0, N) for y odd. 

A chiral ID system violates fermion doublin^^ and 
can onl y be rea lized as an anomalous edge of a ga pped 
2D bnl M^^ I ^^ ^^. The coupled Majorana wire model, \2.A\ 
or figure must therefore also be holographic and living 
on the surface of a 3D bulk superconductor. This can 
be modeled by a stack of alternating layers of spinless 
Px supercond ucto rs (see figure j^a)). The interwire 
backscattering in (2.4) can be generated by bulk inter¬ 
layer electron tunneling and pairing that are not com- 
petin g with the intralayer p ip pairing. Time reversal 
(2.2) extends to the three dimensional bulk by relating 
fermions on adjacent layers. The coupled Majorana wire 
model can also live on the surface of a 3D class Dili topo¬ 
logical superconductor where each chiral Majorana mode 
is bound between adjacent domains with opposite time 
reveral breaking phases 0 = ±7r/2 (see figure 
The discrete translation order along the ^-axis perpen¬ 
dicular to the wire direction can be melted by prolif¬ 
erating dislocations (see figure |^c)). With continuous 
translation symmetry restored, time reversal symmetry 
becomes lo cal w ith = — 1 and the coupled Majo rana 
wire model (2.4) recovers the surface Majorana cone (2.1) 
in the continuum limit for small ky. 



FIG. 3. Coupled Majorana wire model on the surface of (a) a 
stack of alternating px ± ipy superconductors, and (b) a class 
Dill topological superconductor (TSC) with alternating TR 
breaking surface domains, (c) A dislocation. 


The non-local time reversal symmetry (2.2) actually 
provides a weaker topological protection to gapless sur¬ 
face Majorana’s than a conventional local one. For in¬ 
stance in section IIIIl we will show that the N = 2 
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coupled Majorana wire model can be gapped by single¬ 
body backscattering terms without breaking time re¬ 
versal, leaving behind a surface with trivial topological 
order. This reduced robustness stems from the half¬ 
translation component in the ant iferr or magnetic time re¬ 
versal. In the BdG description (2.5), the time reversal 


operator takes the momentum dependent form 


Tk = 


1 


nikn, 


-Ty+i- 


.1 - e 


ik^ 


JC 


( 2 . 6 ) 


normal ordered. The coupled Majorana wire model car¬ 
ries currents that propagate in alternating directions (see 
figure so that Jy{z) are holomorphic for even y and 
are anti-holomorphic for odd y. Focusing on an 
even wire, from the operator product expansion (OPE) 

^ab 

= -+ ... ( 2 . 10 ) 

z — w 

the so{N)i currents obey the product expansion 


for 1C the complex conjugation operator. It commutes 
with the BdG Hamiltonian Tki^g^Q(k) = i7g^Q(—k)Tk 
as well as the particle-hole (PH) CTk = F_kC, for 
C = TxlC the PH operator. In the continuum limit or for 
small ky^T TylC agrees with the conventional local time 
reversal operator and protects a zero energy Majorana 
Kramers’ doublet. The BdG Hamiltonian has a chiral 
symmetry nkifBdG(k) = -iJgjG(k)nk, for Ilk = *CTk 
the chiral operator. It can be used to assign the chirality 
of a Majorana cone by an integral winding number 

n = T / Tr [/i(k)-Wk/i(k)l ■ dl (2.7) 
27rz /c.(ko) 

locally around a loop C£(ko) 5 away from the zero mode 
at ko- Here h(k) is the elliptic operator 


J^{z)r{w) = + W + ... (2.11) 

^ ^ ^ ^ (z-w)^ ^ z-w ^ ^ ' 

where are the structure constants of the soJN) Lie 
algebra with [t^, r] = fp'tsC (see appendix [M. The 
Sugawara energy momentum tensor (along a single wire) 
is equivalent to the free fermion one^ 

T{z) = (2.12) 

for J = (J^) the current vector and -0 = ... ,'0^) 

the A/'-component real fermion. The energy momentum 
tensor defines a chiral Virasoro algebra and characterizes 
a chiral GET. It satisfies the OPE 


h{k) = P+H°^G{k)P^ (2.8) 

for = {P^)‘^ the two local projectors diagonalizing 
the chiral operator Hk = — P^). However, 

as time reversal squares to TkT_k = —which is the 
eigenvalue of the primitive translation —ty at momentum 
k, so does the non-symmorphic chiral operator H^ = 
^-iky ^ The two chiral branches Hk = switch 

across the Brillouin zone when ky ^ kyP27r. As a result, 
a global winding number can only be defined modulo 2. 


A. The so{N)i current algebra 


We notice the couple Majorana wire model (2.4) has 


a SO{N) symmetry that rotates the Wcomponent Ma¬ 
jorana fermion ^y O^'ipy. Gonsequently, there is a 
chiral so{N) Wess-Zumino-Witten (WZW) theor 3 J^^^ or 
affine Kac-Moody algebra at level 1 along each wire. Here 
we review some relevant features of the so{N)i algebra, 
which are well-known and can be found in standard texts 
on conformal field theory (GET) such as Refi53l 

The so{N)i currents have the free field representation 


J^{z) = U{zftP^lz{z) = W r{z)C^\z) (2.9) 


ah 


T{z)T{w) 


c-12 2T{w) dy,T{w) 

(z — w)^ (z — w)^ z — w 

(2.13) 


where the chiral central charge c_ = A/'/2, loosely 
speaking, counts the conformal degrees of freedom on 
the Majorana wires and is proportional to the energy 
currenfP^H^ and entanglement entrop^l^^M^ carried by 
the wire. 

Excitations of the Wcomponent Majorana wire trans¬ 
form acording to the SO{N) symmetry. They decompose 
into primary fields and their corresponding descendants. 
A primary field Va = (E^, •. •, V^) is a simple excitation 
sector that irreducibly represents the so{N)i Kac-Moody 
algebra. 




iti)rs 


S = 1 


w 


V%w)p... 


(2.14) 


where A labels some d-dimensional irreducible represen¬ 
tation of so{N) and is the d x d matrix representing 
the generator of so{N). Eor example it is straightfor¬ 
ward to check by using the definition (2.9) and the OPE 
(2.10) that the Majorana fermion -0 = ..., is 

primary with respect to the fundamental representation. 


i.e. 


where the t^’s are antisymmetric N x N matrices that 
generate the so{N) Lie algebra (see appendix [A ), 2 ) = 
gT+zai complex space-time parameter, and (2.9) is 


6=1 ^ ^ 


(2.15) 
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From ( 2.12| ), space-time translation of a primary field Va 
is governed by 


rrf \^.T f \ \ \{w) 

T{z)Vx{w) = - -^Va(w) H-h . 


{z — wY 


z — w 


(2.16) 


where the conformal (scaling) dimension is given by 

Qa 


hx = 


2{N-1) 


(2.17) 


for — = Qx'^dxd the quadratic Casimir operator. 

For instance the quadratic Casimir eigenvalue for the 
fundamental representation, is — 1 (see appendix [A| ) 
and therefore the fermion -0 has conformal dimension 
= 1/2. This agrees with the OPE ( |2.1Q| ) by dimension 
analysis. 

There are extra primary fields other than the the trivial 
vacuum 1 and the fermion ip- The spinor representations 
(see appendix [A|) a, for N odd, or and 5 _, for N even, 
also correspond to primary fields of so{N)i. Their con¬ 
formal dimensions can be read off from their quadratic 
Casimir values (A7), and are 


^"“16’ ^*^“16- 


(2.18) 


Unlike the infinite number of irreducible representations 
of a Lie algebra, the extended affine so{N)i algebra only 
has a truncated set of primary fields { 1 , a, t/}, for odd, 
or {1, s_, t/}, for N even. 

These so{N)i primary fields take more explicit oper¬ 
ator forms after bosonization and can be found in ap¬ 
pendix and 


B. Bosonizing even Majorana cones 

In the case when N = 2r is even, the N Majorana 
(real) fermions on each wire can be p aired into r Dirac 
(complex) fermions and bosonize^^^ 


Cy 


V2 




exp 


where 0 ^,..., ar e real bosons on the wire, and the 
vertex operator in (2.19) is normal ordered. The bosons 


obey the equal-time commutation relation 


(p{{x),4P,{x') =i7r(-l) 




SyylS^^ sgn(x' — x) 


+ <5!/!/'Sgn(j - j') + sgn( 2 / - y') 

( 2 . 20 ) 

where sgn(s) = s/\s\ = ±1 for s 7 ^ 0 and sgn(O) = 0. 
The first line of (2.20) is equivalent to the commutation 


relation between conjugate fields 

^y{x),dx'^y,{x') =2m{-lY5yy'5^^'5{x - x') (2.21) 


and is set by the “pg” term of the Lagrangian density 




( 2 . 22 ) 


y=-oo j=i 


The second line of (2.20) guarantees the correct anti¬ 


commutation relations between Dirac fermions along dis¬ 
tinct channels. The alternating signs (—1)^ in (2.21) and 
( | 2 . 22 ) specify the propagating directions along each wire, 
R (or L) for y even (resp . odd). Eq.( |2.20| is symmetric 
under time reversal ( 2 . 2 ), which sends 


r4r-i = (-ir4\ 


T(t)lT ^ = (1)1+1 + Try. (2.23) 


We notice time reversal, in this convention, flips the 
fermion parity as it interchanges between the creation 
and annihilation operators. 


The entire coupled Majorana wire Hamiltonian (2.4), 


when N = 2r is even, can be turned into a model of 
coupled boson wires. The total Lagrangian density is a 
combination 


c = Co-n = Co- (^11 + nx) 


(2.24) 


where the Hamiltonian density R = 1-L\\ -\-l-Li_ consists of 
the sliding Luttinger liquicP^^^ (SLL) component along 
each wire 

oo r 

^ll=^x ^ (2.25) 

y= — oo j=l 

and the backscattering component between wires 

oo r 

Hx = -Vy E(-1)""o®(2<+i/2) (2-26) 

y= — oo j=l 

K+i/2=^-^y+i- (2-27) 


The SLL Hamiltonian (2.25) contains the (normal or¬ 


dered) kinetic term i'lpydx'ipy = i{c^ydxCy -h CydxC^y) in 
( |2.4| ) as well as possible forward scattering terms like 
(2 19) density-density coupling {c\,C y){c\,C yY The inter¬ 


wire backscattering Hamiltonian ( |2.26[ ) is identical to 
the second term i'lpy'lpy = i(cjc^+i + c^c^+i) in (2.4). 


This can b e derived directly by applying the bosoniza¬ 
tion (2.19) and the Baker-Campbell-Hausdorff formula 

^ The alternating 


sign (— 1 )^ in (2.26) is crucial to preserve time rever¬ 
sal symmetry (2.23), which relates = 


2^'^+3/2 - 


The r sine-Gordon terms in (2.26) between the same 


pair of adjacent wires mutually commute 




{x') 


= 0 


(2.28) 


and share simultaneous eigenvalues. If there was a sin¬ 
gle pair of counter-propagating wires, these potentials 










































7 


would have pinned + v)^ between 

the two wires. However, they compete with the sine- 
Gordon terms between the next pair of wires due to the 
non-commuting relation 


2^y+l/2(^)>2l?'y+3/2(^') 




(2.29) 


where the unit step function 0{s) = 0 when s < 0, or 
1 when 5 > 0. In other words, the vertex operators 

produces fluctuations to adjacent pairs. 


=2di^^,^{x') + 27r(-l)^0(a;' - x). (2.30) 


The uniform b acks cattering strength Vy^ as protected by 
time reversal (2.2), exactly balances the competing po¬ 
tentials so that the Hamiltonian H = Hy + Hx remains 
gapless. 


III. GAPPING SURFACE MAJORANA CONES 


The previous section describes the gapless surface Ma¬ 
jor ana fermions of a 3D to polo gical superconductor us¬ 
ing a coupled wire model \2A\ . It consists of an array 
of chiral wires, each of which carries N flavors of Ma¬ 
jor ana fermions co-propagating in alternating directions 
(see figure]^. Together with uniform backscattering in¬ 
teractions between adjacent wires, the model captures 
N surface Majorana cones with linear energy dispersion 
about zero energy and momentum (see figure [^. In this 
section we construct explicit fermion interactions that in¬ 
troduce an excitation energy gap to the surface Majorana 
cones while preserving time reversal symmetry. Gener- 
ically, this leaves behind a fermionic surface topological 
order, which will not be discussed until the next section. 

We begin with the simplest case when there are N = 2 
chiral Majorana channels along each wire and correspond 
to two surface Majorana cones. As eluded in section [H| 
due to the non-local nature of time reversal, the coupled 
wire model can be gapped by single-body backscattering 
terms without violating the symmetry. Although this 
cannot be applied to a conventional topological super¬ 
conductor with local time reversal, this model demon¬ 
strates the idea of fractionalization, which can be gen¬ 
eralized to the many-body interacting case and subse¬ 
quently lead to surface topological order. The Hamilto¬ 
nian 'H = 'Ho T^bc consists of the original model (2.4) 
with two fermion flavors 'ij^y = and the inter¬ 

flavor backscattering 


oo 

V-bc =iu (3-1) 

y=-oc 


which is symmetric under the time reversal ( |2.2| ), T : 

^ The BdG Hamiltonian i7BdG(k) = 

^BdG(^) + ^BdG(^) combination of ( |2.5[ ) and 

^BdG(k) [(1 - COS ky)axTz + (1 + cos ky)(JyTy 

-Aiiky{(jyTz ^crxTy)] (3.2) 

which is symmetric under Tk in ( |2.6| ). The energy spec¬ 
trum depends on the relative strength between the two 
interwire couplings +'0y'0y+i) and impy^py^i 

(see figure]^. When rt = 0, the two Majorana cone co¬ 
incide at zero momentum. A finite u separates the two 
until they have traveled across the Brillouin zone and an¬ 
nihilate each other ai ky = ir when u > 2vy. Once an 
energy gap has opened up, the BdG Hamitonian has a 
unit Ghern invariant 

dkx [ dkyTi (Jk) = 1 (3.3) 



where Tr (J),) = Tt {{dk^u'^\dk^ul) - <19*^0) is 

the Berry curvature constructed from the two occupied 
eigenstates below zero energy of i7BdG(k). The 

coupled Majorana wire model t hus behaves like a chiral 
pHip topological super conduct oi®^. However the single¬ 
body Hamiltonian does not possesses a topological order 
in the sense that it does not support anyonic excitations. 
For instance the ^ —ip Z 2 symmetry is global and 
TT-vortices are not quantum excitations of the model but 
rather introduced as classical extrinsic defects. 


U < 2Vy U = 2Vy U > 2Vy 



FIG. 4. Energy spectrum of the N = 2 coupled Majorana 
wire model with inter-flavor mixing. 


This example relies on a simple decomposition of the 
degrees of freedom along each wire, A = 2 = 1 + 1. 
The two Majorana fermions pJy^pJy are backscattered 
independently to adjacent wires in opposite directions. 
Unlike the intra-flavor couplings 'ivy{'ipy'ipy_^i -h 
inter-flavor terms iupjypjy^i freeze independent degrees 
of freedom and they are not competing with each other. 
It is useful to notice that the decomp ositio n breaks the 
SO{2) I symmetry described in section HA, and as a re¬ 
sult the so{2r)i GFT along each wire splits into a pair of 
chiral Ising GFT’s. 

We can now generalize this idea to all A, but with 
many-body interwire interactions. From now on, unless 
specified otherwise, we turn off all single-body scattering 
terms. For instance, the vertical velocity now vanishes. 









Vy = 0, in the kinetic part 1-Lq of the coupled wire model 
(2.4). First we seek a decomposition of the s o{N) i de¬ 
grees of freedom along each wire (see section |nA) into 
a pair of identical but independent sectors (also see fig¬ 


ure 


0 


so{N)i 2 X ^ 


N 


(3.4) 


an energy gap. This is a direct application of the well- 
studied 0{N) Gross-Neveu problenP^H^ in ID. In the 
discrete limit, this is r elated to the Haldane 0(3) antifer- 
rormagnetic spin chairP^^^, the A ffleck - Kennedy - Lieb 
- Tasaki (AKLT) spin chain^^^^ and the SO{n) Heisen¬ 
berg chairP^J*^. When N is odd, the splitting (|3.4| i s less 
trivial. We will make use of the level-rank duality^^^^ 


where are the Kac-Moody subalgebras 


so{n‘^)i D so{n)n x so{n)n (3.10) 


± _ J so{N/2)i for N even 

^ I so{3)s X so for N odd ^ ^ 

to be discussed below. This fractionalization has to 
be complete in the sense that the Sugawara energy- 
momentum tensor exactly splits into 

Tso(N)i = Tg+ + Tg-. (3.6) 

In particular the central charge divides 

c_ {so{N)i) = 2c_ (^iv) = c_ (gjf) + c- (3.7) 


which comes from the fact that the tensor product 
SO{n) (g) SO{n) is a Lie subgroup in SO{n^). In particu¬ 
lar, we will demonstrate the simplest case in section [THB| 
when n = 3. The division of so(9)i can subsequently be 
generalized to so{N)i for all odd N effectively b y wr iting 
= 9 + 2r. This sets = 50 ( 3)3 x so{r)i in (3.4) and 
the corresponding interwire backscattering interactions 

(p). 


A. Gapping even Majorana cones 


and there are no degrees of freedom left behind. Us¬ 
ing the subalgebra current operators J^-±, which are 
quadratic in t/^’s, we construct the four-fermion backscat¬ 
tering interaction 


%.. = « E j'.js; 


y=-oo 


= - E 


y' = — oo 


t2U-1 

* r^,- 


(3.8) 




j2y' + l 


for u positive, and i?, L labels the propagating directions 
of the currents. This is pictorially presented in figure 
andlH 



FIG. 5. Interwire gapping terms ( |3.8^ (green rectan¬ 
gular boxes) between chiral fractional ^^ sectors 
(resp. (g), ©) in opposite direction. 


We begin with the coupled Majorana wire model ( |2.4[ ) 
(or figure Q with N = 2r chiral fermion channels per 


' wire 


and corresponds to the same number of gapless Majorana 
cones. Similar to the previously shown N = 2 case, the 
gapless modes can be removed using simple single-body 
backscattering terms. We however are interested in find¬ 
ing gapping interactions that would support surface topo¬ 
logical order as well. In sect ion [iT^ and appendices |B1 \C\ 
we described the so{N)i WZW theory, which al ong t he 
wire is generated by chiral current operators (|2.9[) 


jM = {-irM 


(3.11) 


We take the alternating sign convention (—1)^ so that 
under time reversal, T• We con¬ 
sider two subsets of generators, so{r)i containing 
for 1 < a < 6 < r, and so{r)^ containing for 

r-\-l<a<b<2r. As they act on independent fermion 
sectors, the two sets of operators commute or equiva¬ 
lently their operator product expansions (OPE) are triv¬ 
ial up to non-singular ter ms. Moreover the Sugawara 
energy-momentum tensor (2.12) for so{N)i completely 
splits into a sum between 


In this section, we design the fractionalization (3.4) 
of 50 (A) 1 for all N and show that the backscattering 


interactions (3.8) open an excitation energy gap without 
breaking time reversal. In GFT context, (3.4) is also 


known as a conformal embedding^^^^^^^. When N = 2r 
is even, there is an obvious decomposition 


5o(2r)i D so{r)i x 5o(r)]^ 


(3.9) 


where the “-h” sector contains 7 /;^,..., whil e the 
one contains the rest ..., . In section 


review how the J 


so(r)I 


* so(r)f 


we 

interactions contribute 


HI A 


Tsoir)t = T^oir)T = -^ E 

a=l a=r+l 

(3.12) 

This ensures all degrees of freedom in 5o(2r)i are gen¬ 
erated by tensor products between those in the so{r)f 
sectors. Precisely this means any 5o(2r)i primary field 
is a fusion channel of the OPE of certain primary field 
pair in 5o(r)j^ and 5o(r)^. Thus as long as the gapping 
terms independently freeze both sectors, they remove all 
gapless degrees of freedom. 
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The backscattering interactions ( |3.8[ ) couples the 
so{r)i sector on the wire with the 5o(r)^ sector on 
the (^ + 1)^^ one. They can explicitly written as 


T^int — ^ 


y= — oo l<a<b<r 


(3.13) 


Firstly, the interactions are time reversal symmetric as 
(3.13) is unchanged by 'ipy Secondly, it 

breaks the 0(2r) symmetry to 0(r)+ xO{r)~. The sym¬ 
metry breaking can be faciliated by forward scattering 
within wires that renormalizes the v elocities differently 
between the 5o(r)^ sectors. Eq.( |3.13 ) is also a combina¬ 
tion allowed by the chiral 0(r) symmetry 








(3.14) 


The chiral symmetry only allows cross couplings 


between adjacent wires. Instead of ( |3.13[ ), an¬ 


other possibility would be its mirror image with sum¬ 
mands This competes with the origi¬ 

nal, but as long as mirror symmetry is broken and their 
strength is asymmetric, an energy gap will open. In the 
following we will ignore the mirror image by assuming it 
is weaker. 


Next we notice that the four-fermion interaction (3.13) 
is marginally relevant when velocity is uniform. The 
dimensionless coupling strength u follows the renormal¬ 
ization group (RG) flow equation 


= +47r(r — 2)u^ 
dX 


(3.15) 


when length scale renormalizes by I eH. This can 
be verified by applying the RG formula among marginal 
operatorJ^ 


^ = -2nY,Cr9m9n 


(3.16) 


where is the fusion coefficient of the OPE OmOn = 
(jrnnQ^ + ... between operators in the perturbative ac¬ 
tion SS = f drdx^^gmOm- In the current case, the 
fusion coefficient OO = —2(r — 2)0 + ... can be evalu¬ 
ated simply by applying the Wick’s theorem of fermions, 
for O = -Ey,a,6^y“V’yVy+iV’y+i- The plus sign in 


(3.15) shows the interacting strength grows at we ak cou¬ 
pling. To show that the backscattering indeed 

opens up a gap, we first focus on a single coupled pair of 
counter-propagating so{r)i channels (see figure |^. 


1. The 0(r) Gross-Neveu model 

Here we conce ntrate on a particular set of backscat¬ 
tering terms in (3.13) at say an even y. We relabel 


= -0^ and for a = 1,..., r. The in¬ 

teraction between the y^^ and {y H- 1)^^ wire is identical 
to that of the 0{r) Gross-Neveu (GN) modeP^^^^ 


'Hgn = “2 (V’i? ■ 


(3.17) 


where the minus sign is from the fermion exchange statis¬ 
tics This GN model is 

known to have an excitation energy gap for r > 2. 

Eor even r = 2n > 2, the Majorana fermions can be 
paired into Dirac ones and subsequently bosonized (see 

) / \/2 ^ e*' 


section 

j 


3- 4/l 

1,..., n. Using 


= (TTl 




R/L) 


^R/L fQj. 


= 4(4)^ + (4)U E (2®^ 

i=i j=i 


(3.18) 


for 20-^ = ^j? ~ (^Iso see (2.27)) are mutually com¬ 
muting variables, the GN interation (3.17) takes the 
bosonized form 

n 

i=i ii/i2 ± 

n 

= u dx4>^jidx4>^L ~ (3.19) 

j = l olEA 

wher e 20 = (20 ^,..., 20^^) and a are roots of so{2n) 
(see j®). The first term renormalizes the velocity W 
in ( |2.25| ) as well as the Luttinger parameter. We assume 
W >> u so that the first term can be dropped. The 
remaining sine-Gordon terms are responsible for gapping 
out all low energy degrees of freedom. Eirstly the an¬ 
gle parameters mutually commute and share simultane¬ 
ous eigenvalues. The ground state minimizes the energy 
by uniformly pinning the ground state expectation value 
(GEV) 


(20-^ (x)) = e Z. 


(3.20) 


We notice in passing that the following subset of sine- 
Gordon terms 


—u cos (a/ • 20) = —u 

1=1 
n 

= -E 


cos 


7=1 


LJ=i 


(nfK#) 


(3.21) 


7=1 


using the simple roots cxj in ( |A9[ ), is already enough to 
remove all low energy degrees of freedom. Here Kjj is 
the Gartan matrix (A12) of so(2n) that appears in the 
Lagrangian density 


To = 

ZTT 


(3.22) 
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for K = K ^ {—K) and ^ = (0^, 0^), and <p is related 
to 0 by the basis transformation ( B13| ). For instance, the 
n vector coefficients nj = (ej, ej) in (3.21) form a null 
basis 


uJKuj = 0 


(3.23) 


and guarantee an energ y gap according to Ref l95l The 
remaining GN terms in ( 3.19[ ) are compatible with (3.21) 

as they share the same minima. _ 

There are constraints on the GEV in ( |3.2Q[ ). In 
order to minimize — cos(a • 20) in ( 3.19[ ), (a • 20) must 
be an integer multiple of 27r. This restricts uniform parity 
among so that the sign in the fermion backscattering 
amplitude 




4(a;)4(a:)^ 


(3.24) 


does not depend on fermion flavor j. This is not the 
only non-zero GEV as is not the only primary field in 
so{2n )i. The backscattering of spinor fields Vs^ = 

( B24| ) corresponds to the two GEV’s 


= (e 




/2 


(3.25) 


where s = (^i,..., 5n) for = ±1, and the overall sign 
Ylj Sj is positive for the even spinor field or negative 
for s_. Here the GEV ( |3.25 ) does not depend on the 
choice of s. This is because given e and s' with the same 
overall parity U^j = £•€) and s' • 0 differ by some 

combination of a • 20, which takes expectation value in 
27rZ. 

There are extra constraints between and from 


the f usion rules of the primary fields of so{2n)i (see (B25) 
and (B26)). Eirstly, s± x requires 


rrisj^ = ms_ + 2m^ mod 4Z. 


(3.26) 


distinct ground statej^ ^^ l ^^ i Eor example, the vertex 
operator ^^^(^o) = of an even spinor field 

creates a jump in the GEV ( |3.24| ) 

(V+(a;o)^e*2Q^(="V4(a;o)) = (3.29) 

because of the Baker-Hansdorf f-Cam pbell formula and 
the commutation relation from (2.20) 


2Q^{x),£\ ■ (j)ji{xo)/2 = m{9{xQ - x)-n^- j -1) 

(3.30) 


for 0 the unit step function 0{s) = 0 when s < 0, or 1 
when s > 0, and m'^ = + n — j + 1. In general, the 

primary fields and corresponds 

to the domain walls of m. 


^s± • 


y4(a;o)+e*4'0Wv;«(a;o) 

4(:ro)te-±-®(")4(:ro) 


_ g¥(™'a±+ 26 '(®o-a:)) 

(3.31) 


Now we move on to the odd r = 2n + 1 > 1 case. 
Eirst we pair the first 2n Majorana fermions into n Dirac 
ones and bosonize them similar to the previous even r 
case. This leaves a single unpaired Majorana fermion 
'^R/L' Dropping terms that only renormalizes velocities, 
the GN model ( 3.17[ ) takes the partially bosonized form 


'Hgn ^ —u 


E 

CXEAso{2n) 


COS (a • 20) 


— u 


E- 

i=i 


cos 


(20^) 




(3.32) 


Take the highest weights 4 = (!,...,!) and £« = I® 

(1,..., —1) for instance, • 0 = • 0 + 20^^ imples 

ms_^{s^) = ms_ (s^) + 2m^. Lastly the fusion rules 


s± X s± 


1, for n even 
-0, for n odd 


requires the GEV’s to obey 


(-ir 

i-iy 


= 1 


for n even 
for n odd 


(3.27) 


(3.28) 


for similar reasons. 

The GN model therefore has four ground states when 
r = 2n > 2. They are specified by the quantum numbers 
(i) rris^ = 0,1, 2, 3 modulo 4 when n is odd, or (ii) rris^ = 
0, 2 and ms_ = 0, 2 mo dulo 4 when n is even. The rest 
are fixed by (3.26) and (3.28). Quasiparticle excitations 


are trapped between domain walls or kinks separating 


and is responsible for gapping out first 2n Majorana 
channels. Projecting onto the lowest energy states and 
taking the GEV (cos(20-^)) = (—1)"^^, the interacting 
Hamiltonian becomes 

T^gn ^ —2n{n — l)u — ni^(—(3.33) 

which is identical to the continuum limit of the quantum 
Ising model with transverse field after a Jordan-Wigner 
transformation. The remaining Majorana channel 
is gapped by the single-body backscattering term. The 
sign of the mass gap nu{ — l)^^ determines the phase 
of the Ising model. We take the convention so that a 
negative (or positive) mass with = 1 (resp. = 0) 
corresponds to the order (resp. disorder) phase. 

Like the previous case, the fermion backscattering am¬ 
plitude (3.24) is n ot the only ground state expectation 
value. Erom (G5) appendix the Ising twist field of 
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so(2n + l)i can be written as the product 

where e = (^i,... ,5^) for Sj = =bl, and 

is the twist field along the last Major ana channel. There 

are three possible GEV for the backscattering 


{Vj^{x)V,^{xy) = (e*-®(-)<7],(:r)<T£(rr)) (3.34) 

J 0 for the disorder phase 
I ±1 for the order phase 


Here we choose the convention so that cfRaL takes the 
role of the spin operator a in the Ising model and its 
non-trivial GEV’s in the order phase specify two ground 
states I t) and | 1). 

Again, quasiparticle excitations are trapped between 
domain walls separating distinct ground statej^ ^^ l ^^ l 
Eor example a twist field (or V^) sits between the 
order to disorder phase boundary where the quantum 
number flips from 1 to 0, or equivalently the fermion 
mass gap in ( |3.33| ) changes sign. This is because the 
twist field V^{xo) introduces a flip in boundary condi¬ 
tion '0i?(xo+) = —iPr{xo—) and corresponds to a change 
of sign in front of the fermion backscattering Al¬ 

ternatively, this can also be understood by identifying 
Va as a Jackiw-Rebbi solitorl^ or a zero energy Ma- 
jorana bound state between a trivial and topological 
superconductoil^ in ID. 

Nex t a t ~ i domain wall of opposite signs of the 
GEV ( |3.34| ) in the order phase traps an excitation in the 
fermion sector This can be seen by equating the order 
Ising phase to a ID topological super conduct oi!^, where 
the two Ising ground states corresponds to the even and 
odd fermion parity states among the pair of boundary 
Majorana zero modes. Adding (or subtracting) a f ermio n 
therefore flips the parity as well as the GEV in (3.34). 
We notice this domain wall interpretation of excitations 
is consistent with the non-Abelian fusion rule 


a X a = 1 (3.35) 

The trivial fusion channel corresponds to the annihilation 
of a domain wall pair such as 


Luttinger parameter. As a result the fractionalization 
(or conformal embedding) so(4)i D so{2)i x 5o(2)i of 
wires with V = 4 Majorana channels does not lead to a 
gapped theory. Instead we turn to an alternative frac¬ 
tionalization so(4)i = 51^(2) X su{ 2 )p^ that only applies 
for V = 4. 

The four Majorana pjy along each wire can be paired 
into Dirac channels Cy = ((ipy -hi'0^)/\/2 = and = 

[tjjy -h ip)^/\f2 = e^^y. It would be more convenient if 
we express the bosons in the new basis using the simple 
roots of so(4): (p^ = (p^ — and = 0^+0^. Unlike 
when r > 2, these bosons decouple in the Lagrangian 
density \2.22\ 

^ oo 2 

^0 = ^ E (-1)" E (3.38) 

y= — oo J =1 

This is equivalent to the fact that the Cartan matrix 
^so( 4 ) = diag(2, 2) is diagonal so that the Lie algebra 
splits into the product su{2)~^ x su{2)~ of isoclinic rota¬ 
tions, each with Gartan matrix Ksu( 2 ) = 2. 

The su{2)i current generators are given by S^{z) = 
i\/2dp^{z) and S^{z) = {S^ ±iSy)/V2 = and 

they satisfy the OPE 

S!{z)Sl{w) = + ... (3.39) 

for / = 1,2 = . The si4(2)j^ sector is completely 

decoupled from the su{2)^ one as the OPE Sl{z)S‘^{w) 
is non-singular. They completely decomposes all low en¬ 
ergy degrees of freedom as the energy momentum tensor 
splits into 

1 2 ^ 

Tso(a), = --Y,d4Piz)d4Piz) (3.40) 

i=i 
2 

= - d^\z)d<t>\z) = 

J=1 


order disorder order 


fusion 


... tt • • •) 


(3.36) 


while the fermion fusion channel corresponds to joining 
the pair of “order - disorder” domain walls into a kink 


order disorder order 


fusion 


H-.-tm---)- (3.37) 


2. The special case: so(4)i = su{2)i x sr^(2)i 


The case when r = 2 requires special attention. The 
0(2) GN model (3.17) is a ga pless Luttinger liquid be¬ 
cause its bosonized form (3.19) contains no sine-Gordon 


terms and the rest only renormalizes velocities and the 


The gapping Hamiltonian is 

oo 

mnt = u Yi ( 3 - 41 ) 

y=-oo 

oo 

= 2u Y 9x4>ldx4>l+i-^cos{4ey+i/2) , 

y=-oo 

4 e,+i /2 = 20^1 - 202 (3,42) 

= + 0y+l +0^-02. 

The first kinetic term of the interacting Hamiltonian only 
renormalizes velocities and the Luttinger parameter. The 
second sine-Gordon term involves four-fermion interac¬ 
tions and is responsible for the energy gap as it back- 
scatters the su{2)p^ sector on the wire to the su{2)^ 
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sector on the {y + 1)^^ one. It pins the ground state 
expectation value (GEV) 




(3.43) 


which characterizes the two distinct ground states. Like 
the previous cases, quasiparticle excitations are kinks in 
the GEV. The fundamental excitation can be created by 
the vertex operator Vs = which is the semionic 

primary field in the su{2)^ sector along the {y + 1 )^^ 
wire. 


B. Gapping odd Majorana cones 

We now move on to the case when there are N = 
2r + 1 > 3 chiral Majorana channels on each wire in 
the coupled Majorana wire model ( |2.4| ) (of figure [^. It 
corresponds to an odd number of Majorana cones on the 
surface of a 3D topological superconductor. The chiral 
degrees of freedom along each wire are described by a 
so{N)i WZW theory, which is goin g to be fractionalized 


into the pair x Q 


iN 


according to (|3.5|). The sector 


along the y^^ wire will then be back-scattered onto the 
sector along the {y + 1)^^ one by the current-current 


interaction (3.8), which will introduce an energy gap. 

Unlike the even N case where so{N)i can simply be 
split into a pair of so(V/2)i’s, here the decomposition 
is less trivial but leads to more exotic surface topolog¬ 
ical order. We begin with the particular case where 9 
Majorana channels can be bipartite into 


so(9)i D so { 3)3 X 50 ( 3)3 


(3.44) 


essentially by noticing that the tensor product SO{3) G 
SO{3) sits inside SO{9). The two 50 ( 3)3 WZW sectors 
carry decoupled current generators. They can the n be 
back-scattered using the current-current interaction (3.8) 


onto adjacent wires in opposite directions (also see fig¬ 
ure and 1^ . 

Eor a general odd V > 9, one can decompose the Ma¬ 
jorana channels into V = 9 + {N — 9). The first 9 chan¬ 
nels can be fractionalized by ( |3.44[ ), which we will dis¬ 
cuss in detail below, and the remaining even number of 
channels can be split using the previous method, namely 
so{N — 9)1 = 50 X so In the case when 

N is smaller than 9, one can add 9 — N number of non- 
chiral Majorana channels to each wire. These additional 
degrees of freedom can be interpreted as surface recon¬ 
struction as they do not violate fermion doublin^^ and 
are not required to live on the boundary of a topological 
bulk. Now each wire consists of 9 right (or left) prop¬ 
agating Majorana channels and 9 — N left (resp. right) 
propagating ones. We still refer the remaining even chan¬ 
nels by so{N — 9)1 except now the negative N — 9 signals 
the reverse propagating direction of these Majorana’s. 

The 50 ( 9)1 and so{N — 9) i sectors can then be biparti- 
tioned independently. The fractionalization of a general 


odd number of Majorana channels is summarized by the 
sequence 

so{N)i D 50 ( 9)1 X so{N — 9)1 ^ X 0^ (3.45) 

for = 50 ( 3)3 X 50 The and ” sectors 

can now be back-scattered independently using (3.8) onto 


adjacent wires in opposite directions. This removes all 
low energy degrees of freedom and opens up an energy 
gap. 


1. The conformal embedding so(9)i 3 so(3)^ x 50 ( 3)3 

As a matrix Lie algebra, 50 (3) is generated by the three 
anti-symmetric matrices I] = (Ex, Ey, Ez) 


Ex = 


/000\ _ /001\ _ /010\ 

00 1 , Ev= ooo , Ez= -ioo . 

Vo -1 oy ’ y V-10 oy ’ V 0 0 oy 


They can be embedded into 5 o( 9 ) by tensoring with II 3 , 
the 3x3 identity matrix, on the left or right 


E+ = E (g) II 3 , E = 1 


E. 


(3.46) 


We denote 5o(3)^ = span{E^, E^, E^} to be the two 
mutually commuting subalgebras in 5 o(9). 

Recall the free field representation (|2.9[) of the 5o(9)i 


WZW current generators ifjV’V2 for an an¬ 

tisymmetric 9x9 matrix, the 5o(3)f current generators 
are given by the substitution of 

J.o(3)±(^) = |r(^)WW (3-47) 

for 2 : = and J = {Jx, Jy, Jz)- Written explicitly, 

J+ = *(t/,23 + ^56 ^ ^89)^ j- ^ -(^47 ^ ^58 ^ ^69) 

J+ = ^ ^46 ^ ^ 79 )^ J- ^ .(^17 ^ ^28 ^ ^ 39 ) 

J+ = J- = + t/, 36 ) 

for Using Wick’s theorem and the OPE 

(w) = — re) + ..., it is straightforward to 

deduce the 50 ( 3)3 WZW current relations 

3(5ii 


xa + —+ • • • (3-48) 

' J [z - wy z-w'^ 


and J^{z)J^{w) is non-singular, for i,j = x,y,z and ^jjk 
the antisymmetric tensor. 


The 50 ( 3) 3 cur rent relations (3.48) differs from the 
50 ( 3)1 ones ( |2.II ) by the coefficient 3 of the most sin¬ 
gular term. This sets the level of the affine Lie algebra. 
The 50 ( 3)3 WZW theory is identical to su{ 2 )q by notic- 
i ng th at the structure factor of su{2) is /jjk = V^^ijk (see 
(3.39) and Rnf l53|) . The su{2) current generators thus 


need to be normalized by 


6^i 


i\/2^ijk 


S^{z)S^{w) = - -L -+ 

' ^ ^ i ^ ^ (^z - w)^ z-w ^ ^ ^ 


(3.49) 
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where the coefficient 6 of the most singular term sets the 
level of the su(2)e affine Lie algebra. 

The Sugawara energy momentum tensors are the nor¬ 
mal ordered product 


^.0(3); 




so(3): 


:(Z). 


(3.50) 


8 so(3)3 

Written explicitly in the fermion representation ( 3.47| ) 
and using the normal ordered product 

(3.51) 

the energy momentum tensor takes the form 

1 9 ^ 

= (^) = -4 E r(z)dr(z) T ^O^(z) (3.52) 

a=l 

0^{z) = 

+ ^4679 ^ ^2356 ^ ^2389 ^ ^5689 ( 3 ^ 53 ) 

for = 'i/;^(z)'ilj^{z)'ip^{z)'ip^{z). The four-fermion 

terms in cancel when combining the “±” sectors, and 
therefore the energy momentum tensor ( 2 . 12 ) completely 
decomposes 


r,o(3)± 


1 ® 

T’-(9). = -^E^“^^“ = ^.o(3)+ + r 


50 ( 3)3 


(3.54) 


a=l 


non-singular, so is the OPE between T rci^+ 


r / Q \ + and (P ^ / Q \ —. 
Each sector carries half the total central charge of 9 Ma¬ 
jor ana channels 


<^so(3)± - 9/4- 


(3.55) 


The primary fields of 50 ( 3)3 = su{2)q are characterized 
by half-integral “angular momenta” s = 0,1/2,..., 3.^ 
Each primary field ..., E/) irre- 

ducibly represents the WZW algebra 

S;{z)vr{w) = j2 (5;X'K“'H + --- (3.56) 


for i = X, y, z and Sf the su{2) generators in the spin-5 
matrix representation. We label the seven primary fields 
by greek letters = 1, <a±, 7±, / 3 , /, each has conformal 
dimension hs = 5(5 + l)/8 (see table|^. In particular 1 = 
Vo is the vacuum and / = V3 is Abelian and fermionic 
with spin 3/2. 

The rest of the primary fields are non-Abelian. They 
obey multi-channel fusion rules 




where the fusion matrix element Vf ^ = 

__ * 1*2 

mined by the Verlinde formulsP^ 

SsiS'Ss2S'Sss' 


■ 

S1S2 / V 




(3.57) 
0,1 is deter- 

(3.58) 


V. 

1 0 + 

7+ 

p 

7- 

a- 

/ 

5 

01/2 

1 

3/2 

2 

5/2 

3 

hs 

0 3/32 

1/4 

15/32 

3/4 

35/32 

3/2 

dg 

1 V 2 +V 2 I + V 2 \/4 + 2^/2 1 + ^/2 V 2 + V 2 1 


TABLE 1. The “angular momenta” s, conformal dimen¬ 
sions hs and quantum dimensions dg of primary fields Vs 
of 50 ( 3)3 = S 1 ^( 2 ) 6 . 


and the modular S'-matri:jP^ 


1 

- sin 
2 


7r(25i + 1)(252 + 1) 

8 


(3.59) 


which is symmetric and orthogonal. Explicitly, the fusion 
rules are given by 


/x/ = l, /X7±=7q=, fxa±=a^, / x ^ =/3 

7± X 7± = 1 + 7 + + 7-, X = 1 -h 7 + 

^x/3 = l-h7++7-+/ (3.60) 

a± X 73- = -h /d, X 73- = + a- -h P 

<^± X /3 = 7 + + 7 _ 

The quantum dimension dg of the primary field V^ is 
defined to be the largest eigenvalue of the fusion matrix 
Ng = It coincides with the modular S matrix 

element iSos/iSoo and respects fusion rules so that 

(3.61) 

S 


They are listed in table [l| 


2. Zq parafermions 


We first study the simplest odd case when there are 
9 Majorana cone s mimicked by the coupled Majorana 
wire model (2.4) with 9 chiral Majorana channels per 


wire. Now that we have bipartite the degrees of freedom 
accor ding to the two 5 o( 3 )f WZW current algebras ir 


(3.47), they can be backscattered independently to ad¬ 


jacent wires in opposite directions (see eq.(3.8) and fig¬ 
ure]^. As the 50 ( 3 ) J sector completely decomposes from 

the 50 ( 3)7 one, the current backscattering /on+ 

^ 50 ( 3)3 M^)I 

between the (^ — 1)^^ and wire does not compete with 
the next pair 

80(3)3 so( 3 )J 

The current-current interaction consists of four- 
fermion terms and is margina lly re levant. This can be 
seen from the RG equation ( |3.16 ) using the operator 
product expansion (J^ • J^+^)^ ^ -t-J^ • (Recall 

the time reversal symmetric convention (3.11) and that 
jyjy rsj i(^—l)yjyp To see that the interaction indeed 
opens up an excitation energy gap, it suffices to focus on 
a single pair of wires with the Hamiltonian 


7/int — 


50 ( 3)3 so(3) 


(3.62) 
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where R/L labels the counter-propagating directions 
along wire y and ^ + 1. 

First we further decompose the so{3)s WZW theory 
by the coset constructiorP^ 


50(3)3 = ^(1)6 X “Ze’’ 




<§0(3)3 _ <§^^( 2)6 

50 ( 2)3 


Therefore, a priori, it is not obvious that this p — a 
decomposition is a splitting of 50 ( 3 ) 3 , and in fact it is 
not. Only the charge sector is entirely belonging 

to 5 o( 3 )^’~ or 50 ( 3 ) 3 ’^. To show this, we go back to 
the energy-momentum tensor ^^ 0 ( 3 )^ in (3.52), say for R 


U{^)6 

(3.63) 


movers. 


where “Ze” refers to the Z^ par afermion CFT model by 
Zamolodchikov and Fatee\P^^. This is done by noticing 
that SO{3) (or equivalently SU{2)) contains the Abelian 
subgroup SO{2) (resp. U{1)) of rotations about the z- 
axis, and on the CFT level, the 50 ( 2)3 WZW sub-theory 
of 50 ( 3)3 (resp. u{1)q C su{2)q) can be bosonized and 
single-out. To do this we first group three pairs of Ma- 
jorana fermions into three Dirac fermions on each chiral 
sector 




(3.69) 


where the total energy-momentum tensor in partially 
bosonized basis is 


' so(9)i^ 








i=i 


+ i’Rdi’R + ^Rd^R + ^Rd^R 


^R — 


'^R + ij’R 

^/2 

i’l + #i 

x/2 


Cp — 


'4’r + i'4’1 

^/2 

i’R + 

^/2 


— 


-h iip' 


(3.70) 


R 


V’l + 

^/2 


and the operator defined in (3.53) is now 


O, 






and bosonize 




^R/L 


Via 


exp 


R/L. 


(3.64) 


-2i 


for j = 1,2,3. The 50 ( 2)3 subalgebra in the R and L sec¬ 
tors are generated by the J~ and currents operators 
in (107) 


350^90^ + ^ ^ ( 

i=i 

cos + cos - . 

cos - (/)5^) t/>: 


(3.71) 




,97 

R 


+ 


= -3i90^, = 3id(l)‘ 


(3.65) 


Eq.(3.71) is dedu ced b y substituting the fermions by 
the boson fields (|3.64|), whose OPE can be found in 
D1|D2|D3 ) n appendix]^ Eor instance, the factor of i in 


(3.71) is a result of mutually non-commuting More 
importantly, 0^, and are completely decoupled. 


where the boson field of the “charge” sector is the average ^s the “charge” sector only appears in T ,( 3 )^.- , it be- 


^R/L 


^\i/L 


^\lL 


^\lL 


(3.66) 


The “neutral” sector is carried by the three boson fields 


longs entirely in 5o(3)^’ . Similarly belongs entirely 
in 50 ( 3 ) 3 ’^. The “Ze” energy-momentum is defined by 
subtracting the decoupled “charge” sector from 50 ( 3 ) 3 . 


_ 






(3.67) 


^so(2)^ ~ ~ 


which are not independent as -h = 0. 

It is straightforward to check that the “charge” and the 
“neutral” sectors completely decouple from each other. 
Eor instance, the Lagrangian density decomposes 


TT. = T 


- - T 


so{3)^’ ^so(2)| 

9 


(3.72) 

(3.73) 


= -\Y.^%drR-lt.9rR^94>‘ 


cr,J 

R 




2^ / .'-^^Rlryt^RlL 

i=i 


(3.68) 


a=7 

1 

2 L 


i=i 


+ 


cos + cos - <P‘^/^ ip' 

cos pp'"/ - (p'^2') t 


I^R 




27r 


i=i 


where the remaining fermions are sup¬ 

pressed, and (—1)^ = 1, (— 1)^ = —1. 

The Lagrangian density (3.68) involves more degrees 
of freedom in 5o(9)f^^ than just 5o(3)^’~ or 50 ( 3 ) 3 ’^. 


and similarly for the L movers. 

The remaining current operators J± = (Jx ± iJy)lV^ 
of 5 o(3)7Jii the R sector and so{3)^ in the L sector 
(see eq.( |3.47| )) now split into “charge” and “netrual” 
parafermion components 


jR/L 




R/L 


(3.74) 
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where the Ze parafermions are given by the combinations the “Ze” and 50 ( 2)3 sectors. 


^ VL + Vl + Vl) (3.75) 

= ^ Vi + e'^^'Vi + Vi) 


^ ^ '^R/L- Unlike the (/)‘^’s, 

here the “neutral” Ze parafermions ^i?/L belongs entirely 
in 5 o( 3 )^’~ or 50 ( 8 ) 3 ’^. This is because and 

both completely sit inside the 50(8)3 ’s as seen above. 
Otherwise one can verified this by comp uting the OPE 
with the energy-momentum tensor ( |8.7Q| ) explicitly 


^ao(3)3^'- 


5/6 

{z — w)‘^ 


^r{w) + 


d^R{w) 

z — w 


+ ... 


(z — re)^ 


-h • . . 

z — w 

(3.76) 




) both 


and r „( 3 )K.+ and T^^^^^R.+ {z)e 

non-singular. Simi lar O PE hold for the L sector. 
The primary fields ( 3.75 ) generate the rest of the Ze 
parafermions (see (D5) in appendix 0 and they obey 
the known Zq structure by Zamolodchikov and Eateev^. 


3. Gapping potential 


Now that we have further decomposed the 5o(3)^ 
currents in each wire into 50 ( 2)3 = U{1)q and Zq 
parafermion components (see eq.(3.74)), the current- 


current backscattering interaction (|3.62|) between a pair 
of wires takes the form of 


"Hint = ^udx(t>%^dx(t>L + 3u 


+/i.c. 


(3.77) 


The first term only renormalizes the velocity of the bo¬ 
son in the 50 ( 2)3 sector. The second term is responsible 
for openning an excitation energy gap. It extracts a Zq 
parafermion T and a quasiparticle from the 5o(3) J 
sector on the wire and backscatter them onto the 
50 ( 3 )^ sector along the {y + 1)^^ wire. This freezes all 
low energy degrees of freedom and the ground state is 
characterized by the Zq expectation value (GEV) 

(^'t(^)^L(^)) ~ = g2^im/6 


for m an integer. 

Like the 0{N) Gross-Neveu model we discussed in sec¬ 
tion IIIIAll quasiparticle excitations here also manifest 
as kinks or domain walls between segments with different 
GEV’s. The primary fields <a±,7±,/5 of the chiral 50 ( 8)3 
WZW theory in table [l| decompose into components in 


a+ = [< 7 i] X a. = [<75] x 

7+ = [ 0 - 2 ] X [e*'^V 7 _ = [ 0 - 4 ] X 

/3 = [ag] X (3.79) 

where ai are primary fields in the chiral Zq parafermion 
theory so that take t he roles of the order param¬ 

eters of the Zq modeP^^. They satisfy the exchange 
relations 

^ix)ai{x') = (Ti(a;')^(a;)e- 2 ” 3 '^U-^') (3.30) 


for R sector, and similar relations hold for the L sector 
with the Zq phases conjugated. Therefore adding the 
operators OL±{x)^'y±{x)^ P{x) to the groun d state create 
kinks of different hights in the GEV (3.78) 


(^al.{xo)'^^ii{x)'^L{x)a±{xo)'j ®o)) 

(7l(^o)4'V^)^i(^)7±(^o)) ~ e^("®±20(.-.o)) 

(y (xo)«'t(a:)^'L(a:)^(^o)) - eTH+ 3 '^U-o)) (3.8I) 

where 0{s) = (5/|5| + l)/2 is the unit step function. 

The fermionic supersector / in 50 ( 8)3 (see tablecon¬ 
sists of operators that admit free field representations. 
Again we focus on the the so{3)^'~ sector. The opera¬ 
tors 


y±2 ^ 

span a 5 = 3 representation of the affine 50 ( 8)3 Lie alge¬ 
bra, where parafermions satisfy¬ 

ing the OPE ^^{z)^^'{w) - (^_^^)-^^73^m+m' 
appendix [P] for explicit definitions). Erom ( |3.8Q| ), they 
create a kink to the order parameter (6) = {0r{x)Pl (x)) 


{Vf{xoy/3R{x)/3L{x)Vf{xo)) = (6)(-l)''U-a.o) ( 3 . 32 ) 


in the order phase. 

The gapping potential can now be generalized to an ar¬ 
bitrary odd number of Majorana channels per wire. Us¬ 
ing the decomposition (3.45), the N Majorana channels 


are first split into 9 + {N — 9). The first 9 channels are 
fractionalized into 5o(3)^ x 5o(3)7 while the remaining 
N — 9 can be split into 5o(^^^^)j^ x 5o(^^^^)7 bec ause 
V —9 is even. The interwire current backscattering (3.8) 
takes the form 


?7int — ^ ^ ^ J 


y=-oo 


y 

■^<^( 3)3 


y?/+l 
so(3) + 


F?/+l 


N-9\- N-9\ + 

soi^ 2 )i ^‘^( 2 )i 

(3.83) 


where different terms act on completely decoupled de¬ 
grees of freedom. They also gap out all low energy de¬ 
grees freedom as the energy-momentum tensor of the 
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CFT along each wire decomposes 


^so(N)i ^so(9)i “1“ ^so(N—9)i 


(3.84) 


- ^.0(3)+ + ^.0(3)3 + ^.o(^)+ + 


using (3.54) and the fact that 

m+n 

2 


^so(m+n)i ^ ^ ^so{m)i ^so{n)i' 

(3.85) 


a=l 


C. Gapping by fractional quantum Hall stripes 



FIG. 6. Gapping N surface Majorana cones by inserting (2 + 
1)D Gn stripe state and removing edge modes by current- 
current backscattering interaction. 


Previously, we designed interwire interactions that gap 
all Majorana modes without breaking time reversal sym¬ 
metry. Here we provide an alternative where each chi¬ 
ral Majorana wire is gapped by backscattering onto the 
edges of two topological stripes sandwiching the wire (see 
figure [^. The topological stripes could be fractional 
quantum Hall states for instance. Similar construction 
has been proposed to describe surface states of topologi¬ 
cal insulatorP^^ 

First we consider inserting between each pairs of Ma¬ 
jorana wire a (2 + 1)D topological state. It supports 
chiral boundary modes which move in a reverse direc¬ 
tion to its neighboring Majorana wire. As adjacent wires 
have opposite propagation directions, the chiralities of 
the topological states also alternates. This alternating 
topological stripe state can be regarded as a surface re¬ 
construction of the 3D topological superconductor. It 
preserves the antiferrormagnetic time reversal symmetry 
(2.2), which relates adjacent topological stripes by re¬ 
versing their chirality. Unlike the coupled Majorana wire 


mode, the topological stripe state itself is a pure (2 + 1)D 
time reversal symmetric system and is not supported by 
a (3 + 1)D bulk. It has a gapless energy spectrum that 
is identical to N surface Majorana cones and is carried 
by the interface modes between stripes (see figure [^b)). 
However the topological stripe state also carry non-trivial 
anyonic excitations between wires. This distinguishes it 
from the coupled Majorana wire model and allows it to 
exist non-holographically in a pure (2 + 1)D setting. 

The Majorana modes along the chiral wires then can 
be backscattered onto the boundaries or interfaces of the 
topological stripes by current-current couplings. In or¬ 
der for the boundary or interface modes to exactly can¬ 
cel the Majorana modes along each wire, the topological 
stripes must carry specific to polog ical orders. We take a 
Gn topological state (see eq.(1.3)) so that its boundary 
carries slQn Kac-Moody current, for Qn the affine Lie al¬ 
gebra of GAT defined in ( |3.5| ) . G^ and G^ denote stripes 
with opposite chiralities. The (2 + 1)D Gn topological 
state itself can be constructed using a coupled wire con¬ 
struction similar to that in Ref l38l and EH and will not be 
discussed here. 

There are two ways the Majorana modes can be 
backscattered onto the topological stripes. The first is 
shown in figure [^a). The N Majorana modes along 
each chiral wire is bipartite into a pair of WZW theo¬ 
ries X according to (3.4). Each WZW theory is 


identical to the CFT along the boundary of an neigh¬ 
boring topological stripe but propagates in an opposite 
direction. It can be then be gapped out by the current- 
current backscattering 


^int = . J 


stripe 

Gn 


(3.86) 


Alternatively, one could first glue the topological 
stripes together (see figure [^b)) so that the line interface 
sandwiched between adjacent G^ and G^ states hosts a 
chiral so{N)i CFT. The stripes can then be put on top of 
the Majorana wire array so that each interface is sitting 
on top of a wire with opposite chirality. The current- 
current backscattering 

nj _ _,Twire Tinterface fo 07^ 

Hint — ^*^soIn)i ' *^so(N)i [O.Q() 

between each Majorana wire and stripe interface gaps 
out all low energy degrees of freedom. 


IV. SURFACE TOPOLOGICAL ORDER 

In the previous section, we described how a coupled 
Majorana wire model, which mimics the surface Ma¬ 
jorana modes of a 3D bulk topological superconduc¬ 
tor (TSC), can be gapped by interwire current-current 
backscattering interaction without breaking time rever¬ 
sal (TR) symmetry. In this section, we pay more atten¬ 
tion to the topological order and the anyon type^^^^^ of 
gapped excitations. The ground states are time reversal 
symmetric and there are no non-vanishing order param¬ 
eters that breaks time reversal spontaneously. There is 
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a finite ground state degeneracy that does not depend 
on s ystem s ize. This signifies a non-trivial topological 
ordepnHinil. 


TR breaking 
gapped surfaee 


so(W)f 




TR symmetrie 
topologieal 
gapped surfaee 




^SO{N)l ■ jfo(A0l or 


— 

- 



- 

L, 1 

3 


Ny 




N 

3^k- 


^_ 

_ W _ 



3 1 

r 

■N_ 

- \ - 




a 


a 


-^ 

< 

FIG. 7. Chiral interface (highlighted line) between a time 
reversal breaking gapped region and a TR symmetric topo¬ 
logically ordered gapped region. 



FIG. 8. The Gn topological order of a quasi-2D slab with 
time reversal symmetric gapped top surface and time reversal 
breaking gapped bottom surface 

The surface topological or der c an be inferred from 
bulk-boundary correspondencd^^^^^^f^. There is a one- 
to-one correspondence between the primary fields of the 
CFT along the (1 + 1)D gapless boundary and the anyon 
types in the (2 + 1)D gapped topological bulk. The con¬ 
formal scaling dimension or spin h = Hr — of a pri¬ 
mary field corresponds to the exchange statistical phase 
0 = corresponding anyon. The fusion rules 

of primary fields are identical to that of the anyons. And 
the modular ^'-matrix of the CFT at the boundary equals 
the braiding 5'-matri:jP^ 

= (4.1) 

in the bulk, where the non-negative integers are the 
degeneracies of the fusion rules 

axb = ^Ar,V (4.2) 

c 

between anyons, and the total quantum di mens ion V = 
A/Sa quantifies topological entanglement^^ and can 


be evaluated by knowing the quantum dimensions da > 1 
of each anyon a by solving the fusion identities 

rfarfb = V^bC^c (4.3) 

c 

On the surface of a topological superconductor, where 
there are no boundaries, the (2 -f 1)D topological order 
corresponds to a (H-l)D interface that separate the time 
reversal symmetric topologically ordered domain and a 
time reversal breaking domain. This interface hosts a 
chiral gapless modes (see figure]^. This geometry can 
be wrapped onto the surface of a slab where the TR sym¬ 
metric and breaking domains occupy the top and bottom 
surface of a 3D bulk (see figure |^. The quasi-2D system 
has an energy gap except along its boundary which is pre¬ 
viously the interface that carries the Qn WZW CFT. The 
bulk-boundary correspondence then determines a bulk 
Gn topological order on the quasi-2D slab. 

Wires in the trivial TR-breaking domain are gapped 
by non-uniform current backscattering 

^TR-breaking = ^ sl{N)i ' ^sl{N)i 

y 

or single-body fermion backscattering perturbation 


^TR—breaking — 1 


to the coupled Majorana wire model (2.4), for A > d 


and ..., i^y). This viola tes the antiferrormag- 

netic time reversal symmetry \2.2\ and leads to a gapped 
surface with trivial topological order. This TR breaking 
half-plane is put side by side against a TR symmetric 
gapped half-plane, where the N Majorana channels per 
wire is fractionalized into so{N)i D x for Qn pre¬ 
viously defined in (3.5). Each Qn sector is then paired 


with the adjacent one on the next wire and are gapped 
by current-current backscattering J^- • Jg+. The inter¬ 
face between the TR-symmetric and TR-breaking regions 
leaves behind one single unpaired fractional Qn chan¬ 
nel. This can be regarded as a 2D analogue of the frac- 
tion al bou ndary modes in the the Hald ane integral spin 
chaiii^IlMI and the AKLT spin chairP^. 

As eluded in the introduction, when the coupled wire 
model involves only current-current backscattering inter¬ 
action, it is a boson model where the bosonic current 
operators, rather than Majorana fermions, are treated as 
fundamental local objects. It is therefore more na tura l 
for us to use the current backscattering Hamiltonian (4.4) 


instead of the fermionic single-body one (4.5) to intro¬ 
duce a time reversal breaking gap. In this case, Tr-fluxes 
are deconfined anyonic excitations realized as 7r-kinks 
along a stripe where there is no energy cost in separating 
a flux-antiflux pair. If the fermionic TR-breaking Hamil¬ 
tonian (4.5) were used instead, 7r-fluxes would be con¬ 


fined on the bottom layer and Majorana fermions would 
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become local. We however will mostly be focusing on the 
former bosonic case, although the fermionic scenario may 
be more realistic in a superconducting medium. 

The bulk-interface correspondence depends on th e ori - 
entation of the time reversal breaking order. In eq.(4.4), 
if the backscattering tunneling strengths are reversed so 
that S > A, figurewill need to be shifted by ^ ^ ^ 1 

and all propagating directions will need to be inverted. 
As a result, the interface CFT will also be reversed to its 
time reversal partner Qjy Qn. This will flip the spins 
of all primary fields ^ = —ha and conjugates all 

exchange phases 0^ ^ 0^ = 0^^. 

An interface with a particular orientation therefore 
corresponds to a time reversal breaking topological or¬ 
der. This is also apparent in the slab geometry in figure 
where the TR breaking order on the bottom surface can 
have opposite orientations. Unlike the conventional case 
on the surface of a topological superconductor where time 
reversal is local, here time reversal involves a half trans¬ 
lation y ^ y -\-l and relates a stripe gapped by J“ * 
to its neighbor • Jy+ 2 - anyonic excitations are 
realized as kinks or domain walls that separate distinct 
ground states along a stripe, time reversal non-locally 
translates anyons on an even stripe (green) to an odd 
one (red) or vice versa (see figure [^. However an inter¬ 
face with a particular orientation can only correspond to 
anyons on stripes with a particular parity. For example 
the bulk-interface correspondence in figure singles out 
anyons on even stripes gapped by J^y ' J%+i- There is 
therefore no reason to expect the anyon theory would be 
closed under time reversal. 


A. Summary of anyon contents 



r even 

r odd 

X 

1 -0 S+ S- 

1 '0 a 

dx 

111 1 

11 

0X 

1 —1 g7rir/8 gTTzr/S 

1 -1 e’""’’''® 


TABLE II. The exchange phase quantum 

dimensions of anyons x in a (2-hl)D SO{r)i topological phase. 

The interface carries chiral gapless degrees of freedom, 
which are captured by the Qn WZW theory whose pri¬ 
mary fields corresponds to the anyon content of the TR 
symmetry gapped surface. For even N = 2r, the surface 
carries a 


Gn = SO{r)i (4.6) 

topological order summarized in table [Ilj Its anyonic ex¬ 
citations obey the abelian fusion rules 

^ X = 1, s± X = Szp (4.7) 

I 1, for r = 0 mod 4 

S± X s± = < 

\ 'ip^ for r = 2 mod 4 


for r even, or the Ising fusion rules 

tp X %p = tp X (j = (j^ (j X (j = 1 -\- %p (4.8) 


for r odd. Eq.(4.7) and (4.8) follows directly from the 


fusion properties of the prim ary fie lds in the so{r)i Kac- 
Moody algebra (see section |II A| and appendix [B| and 
The exchange phase (also known as topological spin) 
can be read off from the conformal dimension 
hx of the primary field Vk in so{r)i that corresponds to 
the anyon type x. Again we exte nd r to negative integers 
by defining SO{—r)i = SO{r)i to be the time reversal 
conjugate of the SO{r)i topological state. 


X 1 


7+ 


/3 


7- 


/ 


1 \/2T72 l + ^/2 yVTT/l 1 + V 2 V 2 + V 2 1 


1 


t /2 ni 


15 + 2r 


X 1 


7+ 


e 16 e 
r even 

p 


— i'Kl‘1 ni 


e 16 


7- 


-1 


/ 


1 V2TV2 I + V2 \/4 + 2s/2 I + V2 V2 + V2 1 




6>x 1 e" 


e 16 
r odd 


TABLE III. The exchange phase dx = ^ and quantum 

dimensions of anyons x in a (2-h 1)D S'0(3)3 ^6 8'0(r)i topo¬ 
logical phase. 


For odd = 9 + 2r, the Gn WZW CFT at the inter¬ 
face corresponds the TR symmetric gapped surface that 
carries a topological order given by the relative tensor 
product 


Gw = ^0(3)3 ^6 ^ 0 (r)i 


(4.9) 


where the fermion pair b = 'ipso( 3)3 ^ '^so(r)i is con¬ 
densed. The concept of anyon condensation^ will be 


demonstrated more explicitly later in section IV B The 
topological state carries seven anyon types and are sum¬ 
marized in table mi For instance, the anyon structure 
matches the primary field content of the so{3)s WZW 
theory (see table |I| when r = 0. The quasipar ticle fusion 
rules of Gat are similar to the 50 ( 3)3 ones in ( |3.6Q[ ) 

/x/ = l, /x7±=7=f, fxa±=a^, / x/3 =/3 

7± X 7± = 1 -h 7 + + 7_, X = 7 + -h 7- (4.10) 

f3 X f3 = 1 7-^ H- 7— + /, X 73 - = -h ol— -h P 


except the following modifications that dependent on r = 
(V-9)/2. 


a± X a± 


1 + 7+, 
/ + 7+: 
/ + 7-, 
1 +7-, 


for r = 
for r = 
for r = 
for r = 


0 mod 4 

1 mod 4 

2 mod 4 

3 mod 4 


a± X 7± = 


+ /d, for r even 
a- + /3, for r odd 


(4.11) 
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This quasiparticle spin and fusion structure will be shown 
later in section [IVB[ The braiding ^'-matrices of the Gn 
states are summarized in appendix [E| 

The Gn se quence ex tends the sixteenfold periodic 
anyon structure^mSHinSl SO{r + 16)i = SO{r)i to a pe¬ 
riodic class of thirty two topological states 


G 


Ar+32 


^G 


AT- 


(4.12) 


This seemingly contradicts the sixteenfold prediction of 
topologically ordered surface states from Ref lTTHTTl This 
is due to the non-local nature of the “antiferromagnetic” 
time reversal symmetry in the coupled Majorana wire 
model. On the other hand, in general there are multiple 
possible gapping potentials that leads to distinct topo¬ 
logical order. For instance, we will show in a subsequent 
section that for N = 16, there is an extended Eg sym¬ 
metry or an alternative conformal embedding that would 
allow a different set of gapping terms but would forbid 
all electronic quasiparticle excitations. 

The thirty two topological states here follow a Z 32 ten¬ 
sor product algebraic structure 


certain quasiparticles that have non-trivially monodromy 
around it .1^ These includes all the tt fluxes 1/2, 3/2 and 
5/2 in the SO{3)s sector, s± (or a) in SO{r)i for r even 
(resp. odd), as well as the tensor product l/2(8)'0, 3/2(8)'0, 
5/2 (g) 7 /, 1 (g) 2 (g) and 3 (g) (or 1 (g) cr, 2 (g) cr and 

3(g)cr). The remaining anyons are local with respect to the 
boson b and survive the condensation, but certain pairs 
are identified if they differ only by the boson condensate, 
a X 6 = a. This includes 3 = 7 /, l(g)'0 = 2, 2(g)'0 = l, 
1/2 0 s± = 5/2 (g) Szf and 3/2 (g) = 3/2 (g) for even 

r, or l/2(g)cr = 5/2(g)cr for r odd. Special care has to 
be taken for the tensor product 3/2 0 cr when r is odd. 
After condensation, the fusion rule of a pair of 3/2 0 a 
becomes 

(3/2(g)(j) X (3/2(g)(j) = (0TlT2T3)(g)(l-|-7/) 

=0+0+l+l+2+2+3+3 

(4.14) 

which has two vacuum fusion channels and indicates that 
3/2 0 (7 cannot be a simple object. This leads to the 
decomposition 


Gati Gn2 — Gni+N2 (4-13) 

where certain maximal set of mutually local bosons from 
Gni and Gn 2 are pair condensed in the relative tensor 
product. We will discuss this in more detail below. 


B. The 32-fold tensor product structure 


We first explain the relative tensor prod uct that de¬ 
fines the Gn topoiogicai state in eq.( |4.9[ ). We begin 
with the tensor product state SO{3)s 0 SO{r)i which 
consists of decoupied SO{ 3)3 = SU{ 2 )q and SO{r)i 
topoiogicai states. The primary fieids of the su{2)q 
WZW CFT are iabeied by seven haif-integrai “spins” 
5 = 0, 1 /2, 1 , 3/2, 2, 5/2, 3 and are summarized in tabie[I| 
and eq.( 3.6Q| ). These correspond to the anyon structure 
of the (2 + 1)D S 0 { 3)3 topoiogicai state. The topoiogi¬ 
cai order of SO{r)i is weii-knowrl^ and was summarized 
earlier in this section. For instance, “spin” 3 corresponds 
to the BdG fermion quasiparticle /, and the half-integral 
“spins” 1/2, 3/2 and 5/2 are Tr-fluxes that give a —1 
monodromy phase of an orbiting fermion. 

In the coupled Majorana wire model where there are 
N = 9+2r Majorana channels per wire, the gapping term 
explicitly seperates the first 9 and finai 2r channeis and 
the current backscattering potentiai does not mix these 
two sectors. This modei wouid therefore give a decou- 
pie SO{3 )3 0 SO{r)i topoiogicai state. However, there 
couid be additionai iocai time reversai symmetric terms, 
such as intrawire forward scattering and 

that mixes the two sectors and condenses the fermion 
pair b = f 30 ( 3)3 ^ ' 0 >so(r)i- f^ct, fermion pair con¬ 
densation is naturai in a superconducting medium where 
the ground state consists of Cooper pairs. The condensa¬ 
tion of the bosonic anyon b resuits in the confinement of 


3/2 0(7 = cr_|_ T cx— 


(4.15) 


where a± are simpie anyons with identicai exchange 
statistics but opposit e ferm ion parity a± x f = and 
obey the fusion ruies (4.11). 



1 

0 + 

7+ 

/? 

7- 

a- 

/ 

r even 

0 

1/2 ® s+ 

1 

3/2 ® s± 

2 

5/2 ® s+ 

3 

r odd 

0 

(3/2® cr)+ 

1 

1/2 ®cr 

2 

(3/2 ®cr)- 

- 3 


TABLE IV. Identification of the seven anyon types in table |IIl| 
as tensor products. 


We summarize the identification of the seven anyon 
types in Gn = SO{ 3)3 SO{r)i as tensor products 
in tabie This explains the exchange statistics and 
quantum dimensions of the quasiparticles in table |m] 


= OJb, 


— dadh 


(4.16) 


with the exception of the non-simple object 3/2 0 a in 
(4.15) where each component a± carries half of its dimen¬ 


sion. The fusion rules in (4.10) and (4.11) are explained 
by the tensor product 

(ai (g) bi) X (a 2 (g) b 2 ) = (ai x a 2 ) (g) (bi x b 2 ) (4.17) 

except in the odd r cases where again the non-simple 
object 3/2 0 a = aj^_\_a- requires special attention. 
The fusion rules (4.11) of a± in the odd r cases are 


fixed by modular invariance. The braiding S'-matrix (4.1) 
is determined by fusion rules and quasiparticle exchange 
statistics. On the other hand fusion rules can also be 
determined by the ^'-matrix using the Verlinde formula 
(3.58). ^Moreover one can define the T-matrix according 


to the quasiparticle exchange statistics 
Tab = ^ab^a 


(4.18) 
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which corresponds to the modular T-transformation in 
the CFT along the boundary. As a consequence they 
satisfies the SL{2; Z) algebraic relatioiPH 

{ST^) ^ (4.19) 

where C- = cr — cl is the chiral central charge of the 
corresponding CFT along the boundary 

c_(Giv) = c_ ( 50 ( 3 ) 3 ) +c_(so(r)i) = ^ + ^ = J- 

(4.20) 


These put a very restrictive constraint on the allo wed 
topological field theory and fix the fusion rules ( |4.11 ) for 
a± when r is odd. The braiding S matrices can be found 
in appendix [E] 

The relative tensor product structure of the sixteen¬ 
fold SO{r)i sequence itself can also be understood using 
anyon condensation 


SO{ri)i ^5 SO{r2)i = SO{n + r2)i 


(4.21) 


where the fermion pair ipi (g) ^2 is condensed. This can 
be verified by a similar condensation procedure as the 
one presented above. For instance, if ri and r 2 are both 
odd, the tensor product ai ( 8 ) (J 2 will become non-simple 
after condensation and decompose into a pair of abelian 
TT-fluxes, 5 + + 5_, with identical exchange statistics but 
opposite fermion par ities s± X 2 p = and are related by 
an anyonic symmetr} F^^FQ 6 [ 

Next we move on to explaining the general relative 
tensor product structure (4.13) of the 32-fold Gn states. 


Eq.(4.21) describes the cases when both Ni and N 2 are 
even, i.e. G 2 ri ^ 2 r 2 — G 2 ri+ 2 r 2 - A similar anyon 
condensation procedure that defined the relative tensor 
product SO{3)s KI 5 SO{r)i above would prove that 


Gn ^6 SO{r)i = GAr+ 2 r 


(4.22) 


for N odd, where the fermion pair b = ^ '0>so(r)i is 

condensed. 

When both Ni = 9+2ri and N 2 = 9+2r2 are odd, each 
of the two Gni = SO{3)s SO{ri)i theories contains 
seven anyon types 1, , 7 ^, /h The tensor product 

state Gni ^ Gn 2 contains three non-trivial bosons 

b = {bo,b+,b-} = {/^ (g)/^,7| (g) 7?.,7i 0 7+} (4-23) 

as 7 ± have conjugate exchange phases 0^^ = ±i. More¬ 
over, these bosons are mutually local. Firstly, bo have 
trivial monodromy around b± as are local with respect 
to the fermion /. Secondly, as there are bosonic fusion 
channels b± x b± = 1 + 6+ + 6_ + ... and b± x bzp = 
6o + b± are local among themselves be¬ 

cause their mutual monodromy phases are trivial. We 
first condensed the Abelian fermion pair bo = ( 8 ) /^. 

The resulting theory contains the following set of (non- 
confined) anyon types 


Gni 


l,/,7i,7^,7i7^,7i7^, \ 

(4.24) 


where some anyon types are identified by the bo conden¬ 
sate, such as / = /^ = /^ and 7 I 7 ?. = 7 + 7 + x bo^ 
and are therefore not listed. Next we condense the non- 
Abelian boson 6 + = 7 + 7 ?., which is already equated with 
= 6 + X 60 . The general condensation procedure of a 
non-Abelian boson was proposed by Bais and Slingerland 
in RefinU In the present case, it begins with the fusion 
theory T of Gwi ^bo Gn 2 f^at only encodes the asso¬ 
ciative fusion content but neglects the braiding structure 
of the anyons. As the boson 6 + is condensed, it decom¬ 
poses as 6 + = 7 + 7 ?. = 1 + ..., which now contains the 
vacuum channel 1. This reduces the fusion theory T 
into a new fusion theory T' ^ where the certain anyons 
in (4.24) become non-simple objects and decompose into 


simpler components while others are identified by the 
boson condensate. This new fusion category T' contains 
the non-confined anyons in the resulting state as well as 
confined non-point-like objects. 


We start with the first line of anyons in (4.24), which 


are all local with respect to the fermion /. The semion 7 ^ 
is self-conjugate as 7 ^ x 7 ^ = 1 + 7 ^ + 7 I. However 7 ?. is 
now also an antiparticle of 7 ^ since 7 + X 7 ?. = 6 + = 1 + ... 
also contains the vacuum channel. The uniqueness of 
antipartner guarantees the identifications 


7+ 






7- 


7- =7+ 


(4.25) 


which obey the usual fusion rules 7 + x 7 + = 1 -h 7 + + 7 - 
and / X 7 + = 7 +. This in turn determines the decompo¬ 
sition of the non-Abelian boson 


6 + = 7 I 7 ?. = 7 + X 7 + = 1 -h 7 + + 7 - (4.26) 


which is consistent with the boson quantum dimension 
dh^ = = 1 + 2d^. Moreover the non-Abelian fermion 

also decomposes 


7 I 7 + = 7 + X 7- = / + 7 + + 7-- 


(4.27) 


Next we move on to the second line of anyons in ( |4.24[ ), 
which are tt fluxes with r espec t t o the fermion /. From 
the ori ginal fu sion rules ( 4.10 ), (4.11) and the identifi¬ 
cation (4.25), (4.26) and (4.27), the tt fluxes satisfy the 
fusion rules 


{ 1 + / + 27 + + 27 _, for ri + r 2 even 

1 + 1 + 7 + + 7 _ + 27 +, for n + r 2 = 3 mod 4 

/ + / + 7+ + 37 -, for n + r 2 = 1 mod 4 

(4.28) 


X = 1 + 1 + / + / + 47 + + 47_ (4.29) 

X {(3^(3^) = 4(1 + / + 27 + + 27 _) (4.30) 

{a\a\) X (a)j_^^) = 1 + / + 87 + + 87 - (4-31) 

{a\al) X (/?i/?2) = 1 + 1 + / + / + + +_ (4,32) 

for A"! = 9 + 2ri and A 2 = 9 + 2r2. 
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These show and must be non-simple be¬ 

cause their corresponding fusion rules contain multiple 
vacuum channels. The decomposition of is simplest 
and applies to all ri, r 2 


+ a\a^_ 


(4.33) 


where a^a‘^ = a\a‘^ x /. For instance, it is straightfor¬ 
ward to check that this decomposition is consistent with 
the fusion rules. and a^l3‘^ are clearly identified 

as they differ only by the Abelian boson bo = f^f^. We 
therefore will simply denote them as a^l3‘^. Moreover, 
one can show that a^l3‘^ and (3^0^ are also identified af¬ 
ter the condensation of the non-Abelian boson 7 ^ 7 ?. = 
1 + 7 + + 7 _ in ( 4.26[ ). This can be verify by equating the 
fusion equations I3‘^) x ( 7 ^ 7 ?.) = x (l+ 7 ++ 7 _). 

The decomposition of a^l3‘^ = depends on the par¬ 
ity of ri + r 2 . 

When ri + r 2 is even, the pair fusion rule for a\a‘^ 
allows it to be simple since there is a unique vacuum 
channel. Moreover as the pair fusion rule is unal¬ 
tered by the addition of a fermion /, it is identical to 
{a\a‘^) X {a\a‘^). This shows Oj^a^ conjugates and 
therefore identifies with which is self-conjugate. 


= a\.a\ = 


In this case, a^/3^ is decomposed into 


1^2 


= cr + a^a 

where we introduce the Ising anyon a that obey 
crxcr = l-h/, a X f = a 

a X 0^0? = 7 + -h 7_, cr X 7± = a^o?. 


(4.34) 


(4.35) 


(4.36) 


The decom posit ion (4.35| ) is consistent with the fusion 
rules (4.31) and (4.29). The reduced fusion category after 
condensing the boson (4.26) is therefore generated by the 


following simple objects 

•^even = (l,7±, 


(4.37) 


when ri+r 2 is even. It has the fusion rules (4.36) together 
with 7 ± X a^o? = cr + . 

When T\ + r 2 is odd, we need to further separate into 
two cases. When r\ +£2 = 3 mod 4, the fusion rule 
of a pair of a^a‘^ in (4.28) forbids it to be simple. It 
decomposes into 

+ 7 + or + 7 _ (4.38) 

where s± are Abelian anyons that satisfy the fusion rules 


S±x s± = l, s±x f 




<§+ X 7 ± = 7 ± (4.39) 


and the fermion parity 7 ± in (4.38) depends on (ri, r 2 ) = 
(0,3) or (1,2) mod 4 but is unimp ortant for the cur¬ 
rent discussion. The deco mposition ( 4.38| ) is consistent 
with the fusion rule (4.28). In this case, the fusion rules 


{a\a\) X in ( 4.31 ) requires a different decompo¬ 

sition of a^l3‘^ than (14.351). 




(4.40) 


The r educed fusion category after condensing the boson 
(4.26) is therefore generated by the following simple ob¬ 


jects 


H = (l>/>s±>7±) 


(4.41) 


when ri r 2 = 3 mod 4. _ 

When ri + r 2 = 1 mod 4, the fusion rule (4.28) 
again forbids a\a\ to be simple. Moreover as the vac¬ 
uum channel is absent, it is no longer self-conjugate but 
instead is conjugate with since it has opposite 

fermion parity and {a\a\) x = l + l + 37 ++ 7 _. 

We decompose 


= 5 + -h 5^+ 


(4.42) 


where s± are Abelian anyons and g± are non-Abelian 
objects that satisfy 


S± X s± = /, s± X / 




5 ±= 7 +xs±. (4.43) 


The decomposition of q;^/ 3^ also needs to be modified 

= 9++ 9-- (4-44) 

One can check that these decompositions are consistent 
with the original fusion rules. The reduced fusion cate¬ 


gory after condensing the boson (4.26) is therefore gen¬ 
erated by the following simple objects 


^'1 = (l,/,s±,7±,fi'±) 


(4.45) 


when ri + r 2 = 1 mod 4. 

Not all objects in the reduced fusion theories 
and 7^3 in (4.37), (4.45) and (4.41) are non-confined 
anyons in the new topological sta tes. S ome may be non¬ 
local with respect to the boson 5+ (4.26) and are therefore 
not point-like objects when 6 + is condensed. They are 
equipped with a physical string or branch cut that ex¬ 
tends. The anyon theory, which encodes both fusion and 
braiding information, after condensation excludes these 
confined extended objects. To determine which objects in 
the reduced fusion categories T' are non-confined anyons, 
we look at the possible monodromy around the condensed 
boson Suppose ai ( 8 )a 2 and bi 0 h 2 a re anyons in the 
tensor product state Gni ^60 ^^2 (4.24) that are related 
by the fusion rule 6 + x (ai (g) a 2 ) = bi (g) b 2 + ..., the 
monodromy under this fixed fusion channel i^^ 


ai (g) a2 



ai (g) a2 



^bi(g)b2 _ ^bi0b2 


Ob+0, 


'ai 


(g)a2 ^ai(g)a2 
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as 6 + is a boson with 6 ^ 5 ^ = 1. In other words triv¬ 
ial monodromy simply reqires the invariance of exchange 
statistics upon an addition of the boson. 

Given any simple object x in the reduced fusion cate¬ 
gory T' in (4.37), (4.45) or (4.41), it may be lifted to mul- 
tiple anyons in the tensor product state Gni ^60 ^N 2 


(4.24) in the sense that it belongs in distinct decomposi¬ 
tions aiG)a 2 = x+... and biG)b 2 = x+.... For instance, 
7 ± are components of the boson 7 ^ 7 ?. = 1 + 7 + + 7 _ 
as w ell as the fermion 7 + 7 + = / + 7 + + 7_ (see (4.26) 
and ( 4.27| )). If x is an object not confined by the bo¬ 
son condensation, then its exchange statistics should be 
independent from the choices of lift 


— ^ai0a2 — ^bi(g)b2 


(4.47) 


since the monodromy (4.46) should be trivial. Other¬ 


wise, the object x has to be non-point-like and extended 
as it does not have well defined statistics. For example 
since j± belongs to the decomposition of a non-Abelian 
boson and fermion, they have to be confined objects after 
condensation. 

The relative tensor product Gni ^6 Gn 2 with the con¬ 


densation of the set of bosons b (4.23) contains non- 


confined anyons in the reduced fusion categories ^even^ 
and 7^3 in (4.37), (4.45) and (4.41). For example 
when ri + r 2 is even, the simple object in (4.37) is 


confined and is not an anyon because it can be lifte d into 
and which have distinct statistics, in (4.35) 


and (4.33). When ri -h r 2 = 1 mod 4, the simple objects 
g± are also confined because they belong in a^l3‘^ and 


which have different spins, in (4.44) and (4.42). 


This shows Gwi is generated by the non-confined 

anyons 


G 


forn+r^even 
1 , for n +r2 odd 


The exchange statistics of a and s± are determined by 
that of their lifts. For instance. 


using table when ri + r 2 is even. This shows 

/ Ni + No^ 

Gn, ^bGN 2 = SO' ' ^ ' 


(4.50) 


when both Ni and N 2 are odd and concludes the 32-fold 
tensor product algebraic structure of the Gw-series. 


V. OTHER POSSIBILITIES 

In the previous sections, we proposed time rever¬ 
sal symmetric interactions that gap the coupled Majo- 
rana wi re m odel and lead to a Gat topological order 
(see eq.(4.6) and (4.9)). The interwire current-current 


backscattering interactions depend on a particular frac- 
tionalization, so{N)i 3 Qn x Qn, of the N Majorana 
channels per wire. However, in special cases, we have 
already seen that alternative decompositions exist and 
correspond to different gapping interactions and topo¬ 
logical orders. For example, at the beginning of sec¬ 
tion m we showed when there are even Majorana chan¬ 
nels per wire, the model could simply b e ga pped by a 
single-body backscattering potential (see ( |3.1[ )) and have 
trivial topological order. This is consistent with the Z 2 
classification of gapless Majorana modes prote cted by the 
“antiferromagnetic” time reversal symmetry ( 2 . 2 ). An¬ 
other example was given in section [III A 2 | for the spe¬ 
cial case when there are A" = 4 Majorana channels per 
wire where the decomposition needs to be changed into 
so{4:)i Z su{2)i X 51^(2) 1 . The resulting gapped state 
carries the SU{2)i semion topological order instead of 
G4 = 5'0(2)i. 

Moreover the sixteenfold classification of topolog¬ 
ical superconductors (TSC) with the presence of 
interact iorP^"^i^ suggests the 32-fold Gw-series could have 
redundancies. On the other hand, the Ziq classifica¬ 
tion of TSC is based on the canonical local time rever¬ 
sal symmetry, which is fundamentally different from the 
non-local “antiferrormagnetic” time reversal considered 
in this manuscript. The Z 32 structure of surface topo¬ 
logical order could be an artifact of such unconventional 
time reversal symmetry. Nonetheless, here in section [V A| 
and VB[ we discuss altenative gapping interactions when 
A = 16 that removes all electronic quasiparticles. 


A. Consequence of the emergent Eg when A = 16 

We design alternative inte rwir e backscattering terms 
in the coupled wire model (2.4) with A = 16 Majo¬ 


rana channels per wire. They open a time reversal sym¬ 
metric energy gap among 16 surface Majorana cones 
with the same chirality. In general, these terms can 
also apply when the number of chiral Majorana chan¬ 
nel per wire is larger than 16 by acting on a subset of 
channels. We begin with the bosonized description pre¬ 
sented previously in section ing where each wire consists 
of an 8-component chiral 1/(1) boson 4> = 
that bosonizes the complex fermions Cj = (' 02 j-i + 
iV’2i)/V2 = exp(i0-^). This theory is special because it 
carries non-trivial bosonic primary fields, which can con¬ 
dense. For example the two spinor representations s± 
correspond to bosonic primary fields of go(16)i with con¬ 


formal dimension hs^ = 1 (see eq.( 2.18| )). In particular 


we will focus on the even sector It consists of vertex 
operators 


S + 




S = (5i,...,58) 


(5.1) 


(see eq.( B24| )) for sj = ±1 with si...ss = +1. The 
128 = 2‘ number of combinations naturally matches with 
the dimension of the even spinor representation of so(16) 
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(see appendix IA|. These are related to each other 
through the OrE with the raising and lowering opera¬ 
tors Qf so(16)i (see ( |B8[ ) in 

appendix [B|). The 128 lattice vectors s/2 extend the 112 
roots OL of so(16) to the root lattice of the exceptional 
simple Lie algebra with size 240The unit dimen¬ 
sional vertex operators themselves can be regarded 
as raising and lowering operators that enlarge the so(16)i 
current algebra to E^ at level 1. This extends the current 
algebra of each wire 


5o(16)i C {Es)i (5.2) 


and is intimately related to the fact that the surface state 
can be gapped out without leaving electronic quasiparti¬ 
cles which are non-local with respect to the boson 

The gapping strategy is to condense primary fields in 
the bosonic sector between adjacent wires. This is 
facilitated by interwire backscattering interactions that 
bipartite the emergent E^ symmetry. 

3 5o(8)j^ X 5 o(8)^ (5.3) 


However, these so{S)i subalgebras are distinct from the 
ones in the decomposition so(16)i 3 so{S)i x so{S)i. In 
particular, we will see that they do not support electronic 
prim ary fields Cj = . Out of 128 s lattice vectors in 

(5.1), there is a (non-unique) maximal set of 8 orthonor¬ 
mal vectors s^i),..., S(8) 


1 

2^{m) 



= 25 


mn- 


(5.4) 


We choose the set containing the highest weight vector 

£(i) = (1,1,1,1,1,1,1,1); 


1 

1 \ 


/ 1 1 1 1 -1 -1 -1 -1 \ 

1 1 

|_ 

11 - 1-11 1 - 1-1 
11 - 1 - 1 - 1-11 1 


£( 1 ) • 

1 

^( 8 ) 


1 - 11 - 11 - 11-1 

1 - 11 - 1 - 11-11 

1 

1 J 


V 1 -1 -1 1 1 - 1 - 11 / 


From (2.20), they give 8 mutually commuting bosons 
£(n) • 0S//2 per wire 


2 ^{m) ■ 4^y (^5 1) ) 2 

= 2'Ki5mn{-'^Y5yy'Sgn{x' - x) (5.6) 


up to a constant integral multiple of 27rL 

We separate the 8 vectors into two groups 5+ = 
{^(1)5 ^(2)5 ^(3)5 ^( 4 )} and S~ = {s( 5 ), S( 6 ), S( 7 ), S( 8 )}. 

They defines the two 50(8)f subalgebras in ^8, whose 
roots lie in the root lattice of E^ orthogonal to 5^ re¬ 
spectively. One could pick the simple roots 

OLi = S(i)/2, 0L2 = E 02, CX2 = 63 + 64, = 65 + 66 

= ^( 5 )/ 2 , OL2 =e2 - ei, OL2 = e^- 63, 54 = 66 - 65 


so that their inner product recover the Cartan matrix of 
50(8) 


nf ■ = Ku, K=(-_l ~i / ) (5.7) 

\-l 0 0 2 / 

while opposite sectors decouple • a j = 0. 

The new gapping potential is constructed by backscat¬ 
tering the two decoupled so(8)f currents to adjacent 
wires in opposite directions. 


CX3 




(5.8) 


However not every terms can be written down as 4- 
fermion interactions. In particular Hint contains inter¬ 
wire 5 + quasiparticle backscattering 


+ h.c. - cos 




(5.9) 


for Sj^e'j = ±1, that condenses pairs of s+’s along ad¬ 
jacent wires and confines all electronic excitations. The 

so(8)^ WZW CFT carries three emergent fermionic pri¬ 
mary fields 




= exp 


1 

2 





(5.10) 


for p = 1, 2,3. All of which have neutral electric charge 
and even fermion parity with respect to the original 


electronic operators Cj = e*' 


This is because the 
c^’s are invariant under the U{1) gauge transformation 
E p. As a result, the interaction (5.8) corre¬ 


sponds to a gapped SO{S)i topological order but con¬ 
tains no electron-like anyon excitations. Lastly we notice 
that this matches with the surf ace top ological order of a 
type-H topological paramagnet! ^^ * ^^^ 


B. Alternative conformal embeddings 


The fractionalization so(9)i ^ so{3)s(S)So{3)s in 
is the corner stone for the construction of symmetric gap¬ 
ping interactions when there is an odd number of Majo- 
rana species. However, this is not the unique decompo¬ 
sition. In general when the number of Majorana chan¬ 
nels is a whole square, the wire can be bipartitioned into 
5o(n^)i 3 so{n)n 0 

For instance, this provides yet another alternative 
when N = 16 where each wire is fractionalized into a 
pair of 50(4)4 = 544(2)4 X 544(2)4. The 5 o( 4 )J current op¬ 
erators can be constructed in a similar fashion as those 
in the so{3)f case, J = for = 5] (g) 14 and 

= 14 (g) 5] where 5] are antisymmetric 4x4 matrices 
generating 5 o( 4 ). After introducing the current-current 


HIBl 
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backscattering interactions , , the surface 

would carry a 5'0(4)4 = SU{ 2)4 x SU{ 2)4 topological 
order. Each SU{2)4 theory contains five anyon types 
j = 0 ,1/2,1, 3/2, 2 with spins hj = j{j + l)/6. The 
SO (4 :)4 topological state does not carry fermionic exci¬ 
tations, and therefore, like the previous example in | V A[ 
this gapping potential also removes all electronic quasi¬ 
particle excitations. 

The gapped symmetric states for N odd are not 
unique either. For example, the decomposition so(25)i ^ 
50 ( 5)5 ^ 50 ( 5)5 leads to a surface 5'0(5)5 topological or¬ 
der which is inequivalent to G 25 = SO{3)s KI 5 SO{S)i. 


VI. CONCLUSION AND DISCUSSION 


We constructed a coupled Majorana wire model in (2 + 
1)D that imitates the massless Majorana modes on the 
surface of a topological superconductor. This model had 
a non-local “antiferromagnetic” time reversal symmetry 
and consequently was Z 2 classified - rather than Z in the 
class Dill TSC case - under the single-body framework. 
Despite the difference, this model adequetly described 
the surface behavior of a TSC when the number N of 
Majorana species was odd, and it was worth studying 
and interesting in and of itself. 

We introduced the 4-fermion gapping potentials in sec¬ 
tion in They relied on the fractionalization or biparti¬ 
tion of the so{N)i c urren t alo ng ea ch wire into a pair of 
Qn channels (see eq.(3.4) and ( |3.5[ )). The two fractional 
channels then were backscattered onto adjacent wires in 
opposite directions. This freezed all low energy degrees 
of freedom and opened an excitations energy gap with¬ 
out breaking time reversal symmetry. When N = 2r 
was even, each wire could simply be split into a pair of 
Qn = 5o(r)i channels. The fractionalization was not 
as obvious when N was odd. We first made use of the 
conformal embedding that decomposed nine Majorana’s 
into two su bsectors, 5o(9)i 3 50 ( 3)3 50 ( 3)3 (see sec¬ 

tion IIIBl). This division could be generalized by all 


odd cases by splitting a subset of 9 Majorana’s into a pair 
of 50 ( 3)3 and the remaining even number of Majorana’s 
into a pair of 5o(r)i. This could even be applied when N 
is less then 9 because each wire could be reconstructed by 
adding an arbitrary number of helical Majorana modes 
with the same number of right and left movers. 

The surface Gn topological ordered was infer red from 
the bulk-boundary correspondence (see eq.(1.3)). These 
topological states followed a 32-fold periodicity Gn — 
Gn -\-32 and a relative tensor product structure Gni ^6 
Gn 2 — Gati+ats- We presented the quasiparticle types as 
well as their fusion and braiding statistics properties. We 
explained the relative tensor product structure using the 
notion of anyon condenstiorP^. On a more fundamental 
level, one should be able to deduce the topological order 
without the knowledge of the boundary by studying the 
modular properties of the degenerate bulk ground states 
under a compact torus geometr}®^, or by directly looking 


at exchange and braiding behaviors of bulk excitations. 
In fact the coupled wire construction provided a fitting 
model for this purpose. Being an exactly solvable model, 
a ground state could be explicitly expressed as entangled 
superposition of tensor product ground states between 
each pair of wires. In the simplest case when the model 
is bosonizable, a ground state could be specified by the 
pinned angle variables of a collection of sine-Gordon po¬ 
tentials. The bulk excitations could be realized as kinks 
between a pair of wires and could be created by vertex 
operators. The virtue of a bulk description is that the 
action of time reversal on quasiparticle excitations could 
be examined explicitly, which we have not performed or 
addressed here. These issues are beyond the scope of 
this article and we refer a more detail discussion to sub¬ 
sequent works. 

We noticed that there were alternative ways of frac- 
tionalization that led to different gapping interactions 
and consequently different topological orders. We saw in 
section |III A2| that N = 4 was an exceptional case that 
requires the special bipartition 5o(4)i Z su{2)i x su{2)i 
instead of two copies of 5o(2)i. We also saw in section M 
that when V = 16, the surface could be gapped by al¬ 
ternative interactions that corresponded to a SO{S)i or 
30 ( 4)4 topological order, none of which contained elec¬ 
tronic quasiparticle excitations. Other conformal embed¬ 
dings 5o(n^)i 3 so(n)n G so(n)n could give rise to mul¬ 
tiple possibilities. Our 32-fold topological states, which 
only utilized 5o(9)i ^ 50 ( 3)3 O 50 ( 3 ) 3 , therefore should 
belong into a wider universal framework. These should 
be addressed in future works. 


Appendix A: The so{N) Lie algebra and its 
represent at ions 


The so(N) Lie algebra are generated by real antisym¬ 
metric matrices with entries 

V J NxN 

= (Ai) 

for r, 5 = 1,..., V. There are N(N — l)/2 linearly inde¬ 
pendent generators since and = 0. In 

the main text, we write the basis labels as /d = (r5), for 
r < 5 , for conciseness. The generators obey the commu¬ 
tator relation 


firs) ^^(pq) 


firs)(pq)(mn)t^ ^ 

m<n 


where the structure constant is 


(A2) 


J (rs)(pq)(mn) ^mr^nq^sp ^mr^np^sq 

ms^rq^np ^ms^nq^rp- (A3) 


The matrix representation (Al) is referred as the fun¬ 
damental representation of so(N) and is labeled by 
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In general the generators of so{N) can have different ir¬ 
reducible matrix representations labeled by A. 

Since the quadratic Casimir operator 


The integer r is the rank of the so{N) Lie algebra and 
is determined hy N = 2r for N even or N = 2r 1 for 
N odd. These generators satisfy 


2a = -E^a^a (A4) 

/3 


commutes with all the generators, it must have a fixed 
eigenvalue Qx that (incompletely) characterizes the ir¬ 
reducible representati on A . For instance, the fundamen¬ 
tal representation in ( |A1| ), denoted by has quadratic 
Casimir value = N — 1. 

The spinor r epre sentation a of so{N) makes use of the 
Clifford algebrsP^ {7a, 76} = 7a76 + 767a = ‘^^ab where 
7i, • • •, 7n are hermitian matrices of dimension d = 2^/^ 
for N even or d = for N odd. The so{N) gen¬ 

erators are represented as the quadratic combination 


ab 


(A5) 


and satisfy (A2). When N is even, the parity operator 


( — 1)^ = .. .yAT commutes with all and the 

representation is decomposable into cr = 5+ 0 s_, where 
are 2^/^“^-dimensional sectors with (—1)^ = ±1. 
The so{N) generators are then irreducibly represented 
by 


[H\E^] = [E^,E-^] 





(AlO) 


cx 


ifa+;3eA 

0, if otherwise 


for OL ^ (3. 


The Cartan matrix K = {Kij)rxr of the algebra is 
defined by the scalar product 


2ol^cxj 2odjodj 


(All) 


so{2r) is simply-laced in the sense that all roots have 
identical length and the Cartan matrix is therefore sym¬ 
metric 


^so{2r) 


/ 2 -1 0 ... 0 \ 

-1 2 : 

0 2 -1 -1 

: -12 0 
\ 0 ... -1 0 2 / 


(A12) 


^ = P±t^EPl (A6) 


where P± are the projection operators onto the fixed par¬ 
ity subspaces. As = ~(f/4)l, the quadratic 

Casimir values ( |A4[ ) of spinor representations are 


Qa 


N{N -1) 

8 


Qs± 


N{N -l) 

8 


(A7) 


The complexified so{N) Lie algebra has an alternative 
set of Cartan-Weyl generators. It consists of a maximal 
set of commuting hermitian generators ..., , and 

a finite set of raising of lowering operators E^ = {E~^y ^ 
labeled by integral vectors a = ..., a^) G A called 

roots. The root lattice is given by the set 


^so( 2 r) = {±6/ ±ej:l</<J<r} 

^so(2r+l) = ^so(2r) U {±6/ : 1 < I < v} (A8) 

where e/ are unit basis vectors of W. In particular, there 
are r simple roots ai,..., that forms a basis for the 
root lattice. For so{N) they can be chosen to be 


Sometimes it would be convenient to use the Chevalley 
basis so that the commuting generators are redefined 


= 


\o^i\ 




(A13) 


i=l 


so that the commutator relations ( |A10 ) becomes 

=±KijE^°‘X =5" hr 


(A14) 


Appendix B: Bosonizing the so(2r)i current algebra 


Here we review the bosonizatiorP^J^^HH of a chiral wire 
with N = 2r Majorana fermions, and express the so{2r)i 
current operators in bosonized form. The 2r Majorana 
(real) fermions can be paired into r Dirac (complex) 
fermions and bosonized into the normal ordered vertex 
operators 


c>{z) 


^ (z)+ 
V2 



(Bl) 


Oil 


e/ - e/+i, for / = 1,.. .,r - 1 
< Gr, for / = r and N odd 

er-i + e^, for / = r and N even 


(A9) 


Here we focus on a single wire, say at an even so that all 
fields depend on the holomorphic parameter 2 : = 

The r-component boson 4> = (0^, • • •, 0^^) is governed by 
the Lagrangian density 


The set of roots A consists of integral combinations of \ \ ^ ^ 

the simple roots a = X}j=i ^bat its length is Co = ^ ^ = -^dx(j>dt4> (B2) 

|a| = \/2, for even A, or |a| = 1 or \/2, for odd N. 
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and follows the algebraic relations 


4r>{x,t),Y {x',t) 

= iir 

pj sgn(x' — x) + sgn(_) — j') 



(B3) 


or equivalently the time-ordered correlation function 


(w)) = -5^^' log( 2 : -w) + ySgnO' - f) (B4) 

for sgn(s) = s/\s\ when 5 7 ^ 0 and sgn(O) = 0. Operator 
product expansions between unordered vertex operators 
can be evaluated by 

for B linear combination of the bosons . For in¬ 
stance, the vertex operators in (Bl) reproduce the prod¬ 


uct expansion of a pair of identical Dirac fermions 


c^{z) (c-^(re))^ 


1 


-h idcj)^ {w) + ... 


(B5) 


and the singular piece is dropped when the product is 
normal ordered in the limit z ^ w. The non-singular 
sign factor i7rsgn(j — j') ensures fermions with distinct 
flavors anticommutes 

{z)c^ (w) = —c^ {w)c^{z). (B 6 ) 


The so{2r)i currents in the Cartan-Weyl basis can now 
be bosonized 


or (ii) a • /3 = — 1 in other words a + /3 G A. To sum¬ 
marize, the Cartan-Weyl generators satisfy the product 
expansion 

HHz)HUw) = , - d^Hw)d^{w) + ... 

[z — wY 

H\z)E°‘{w) = -^E°‘{w) + ... 
z — w 

E^{z)E-^{w) = - —(B12) 
(z- wY ^ z-w 

' ' i=l 

- i (a • d^{wp +... 

E‘^{z)El^{w) = '^AA1 e‘^+I^{w) + ..., if a • = -1. 

z — w 


For instance, the 2-cocyle coefficient e{oL^f5) ensures the 
OPE between E^{z) and E^{w) commute as the sign in 
(Bll) when exchanging ol ^ (5 cancels that in l/{z — w) 
when switching z ^ w. 


In certain derivations, especially when involving quasi¬ 
particle excitations, it may be more convenient to use 
the Chevalley basis. Here fields are expressed in terms of 
non-local bosons cp = ( 0 ^,..., 0 ’^), which are related to 
the original ones by the (non-unimodular) basis transfor¬ 
mation 


HYz) = (A {z)(A (z)^ = idzPYz) (B7) 

r 

E^{z) = = exp (^ioL • cf>{z)^ 

i=i 


= (B13) 

7=1 


wher e ol = ..., G A are roots of so{2r) (see 

(A 8 )) and the fermion products are normal ordered. For 
instance, a has two and only two non-zero entries and 
E^ must be of the form 


E°^{z) = c\z)^Y{z)^ = 

Combining raising or lowering operators give 


(B 8 ) 


using the simple roots otj = {a},... G (see (A9) 
appendix 0- The Lagrangian density ( |B2[ ) now be¬ 


m 


comes 


£0 = T ^ Kijd,(t>^dY 


I,J=^ 


(B14) 


E'^{z)El^{w) 




gi(«-0(2;)+/3-0(w)) 

{z — 


(B9) 


where the vertex operator here is again normal ordered 
and the 2 -cocyle is given by the star product 


where K = {Kj j)rxr = 
so{2r)i (see eq.(A 12 )). 


Oil • OLj is the Cartan matrix of 


The current generators are rewritten in the Chevalley 
basis by 


e(a, (3) = (- 1 )“*^ = . (BIO) 

As Yll=i even for all roots, we have the following 
simplification when interchanging ol ^ (5 

e(a,/3)e(/3,a) = (-1)“''®. (Bll) 

Using the boson OPE ( |B4| ), the product of the two ver¬ 
tex operators above is singuiar only when (i) ol = —/3, 


hi{z) = Y(^\H\z) =iYKijdA^{z) 

i=l J =1 

E'^(z) = E^(z) = exp (^ih'^xYiz)) (B15) 

where jd = ^jb'^cxj are roots expressed in integral com¬ 
binations of the simple ones, for b = ( 6 ^,..., 6 ^^) G 
The Chevalley generators satisfy the modified current re- 
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lations from (B12) 


hi{z)hj{w) = .2 + • • • 

[z — wY 

E'°{w) + 


hi{z)E^{w) = 


E'^{z)E-'^{w) 


z — W 

1 

{z — wY 


(B16) 




+ V- hi{w) + ... 

z — w 


TZ^2fr,n\ — ^ 2 ) ^bi+b2 


E'^^{z)E'^^{w) = 


z — w 


I=l 

E'^^+*'^{w) + ... 


if hjKh2 = -1. 

The (normal ordered) energy-momentum tensor can 
be turned from the Sugawara form (2.12) to the usual 
bosonic one 


so{2r) at level 1. When the exponenet a-b in (B22) is —1, 
the vertex operators and are related by the 

SO{2r)i symmetry and belong to the same primary field 
sector. For example the unit vector a = ei is the highest 
weight that generates the fermion sector Y- Applying 
lowering operators E~^ to gives all 2r Dirac 

fermions 


= span 




j = l,...,r} (B23) 


whic h in t urn irreducibly represent the so{2r)i algebra 
(see ( 2.14[ )) according to the fundamental vector repre¬ 
sentation. 

The unit vectors a = e^_i and generate the two 
spinor sectors S- and 5 + respectively. Each of them con¬ 
sists of 2'^~^ twist fields 


T{z) = 


1 


2(77-1) 


'^H\z)H\z) + E‘^{z)E-°‘{z 

oteA 


.i=l 


= • dcf){z) = - ■ Kd4){z). (B17) 


=a'^ ...a 


2r 


(B24) 


= span exp I i y] 1 ' b 


i=i 


i=i 


Excitations in the GET can be easily represented by ver¬ 
tex operators 

V^(z) = exp (ia • (t>{z)) = exp ^iav • 0(^)^ (B18) 

labeled by integral lattice vectors a = (ai,...,a^), or 
equivalently dual root lattice vectors ay = (a^,..., a^) 
with rational entries 

al = (B19) 

IJ 


The conformal dimension of can be read off by the 
inner product 


h^ = Y'^K-^a=^{K-^Y^aiaj 

~ 2^v^v = (B20) 


This can be evaluated fro m de finition (2.16) using the 
energy-momentum tensor ( |B17| ) and the OPE 


dz(l)i{z)(j)j{w) = -{K ^)^'^log(2: - re) + ... (B21) 


which is equivalent to (|B4 


Most vertex operators 


primary and do not represent the so{2r)i 
algebra. The OPE with the current generators 


B18|) however are not WZW 
Kac-Moody 


hi{z)V^{w) = -^v^{w) + ... 
z — W 

E'={z)V^{w) = cl{z - w)^-'°V^+^'°{w) + ... (B22) 


would match the requirement (2.14) for a primary field 
only when the exponent of the singular term is bounded 
below, i.e. a • b > —1 for all roots (3 = Such 

lattice vectors a are called weights or Dynkin labels of 


They irreducibly represent the so{2r)i algebra according 
to the even and odd spinor representations. These are 
the only primary fields of so{2r)i and their conformal 
dimensions are given by = 1/2 and hs^ = r/8. 

The four primary fields 1,7/, obey a set of fusion 
rules, which are OPE keeping only primary fields. 

s± X Y = (B25) 

1, for r even 
t/, for r odd 

(B26) 



t/, for r even 
1, for r odd 


Eor instance, the OPE 


(X [z — w) 2 + ... 

= {z-w)-Ws_{w)E... (B27) 


shows X 7 / = 5_, and 


eiEj7"W/2g-iEy"W/2 oc + ... (B28) 


shows 5 + X = 1 for r even, or x 5_ = 1 for r odd. 


Appendix C: Bosonizing the so(2r + l)i current 
algebra 

A chiral wire with A = 2r -h 1 Majorana fermions can 
be partially bosonized by g roup ing 7 /^,..., in pairs to 
form r Dirac fermions (see ( |B1[ )). This leaves a single Ma¬ 
jorana 7 /^^+^ behind. In order for the fermions to obey 
the correct anticommutation relations, the bosonized 
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complex fermions (Bl) have to be modified by a Klein 
factor 

where (—1)^ is the fermion parity operator that anticom¬ 
mutes with and both 11 and ' 02 r+i commute with 

the rest of the bosons (j)^. In a non-chiral system, (—1)^ 
can be chosen to be the combination for 7l/r fhe 

zero mode of In the chiral case, it can be defined 

by i77oo using an additional Major ana zero mode 7oo 
that completes the Cliffort algebra { 7 , 700 } = 0. 

The so{2r + l)i current algebra extends th e so {2r)i 
algebra by the short roots with length 1 (see ( |A8[ )). It 
contains the so{2r)i genera^rs = id(l)^ and = 

Aoc-cf) 


(see 


m 


apendix |b|, for ol G fbe long 

roots with length |a| = a/2. The remaining raising and 
lowering operators with the short roots are represented 
by the normal ordered products 




(C2) 


In addition to (B12), the Cartan-Weyl generators satisfy 
the current relations 


W{z)E^^J{w) = - E^^i{w) 


E^i{z)E-^i{w) = 


- W 

1 


+ 


-W{w) (C3) 




{z — wY z — w 

— -50-^ (rc)90-^ (re) 

— {w) + ... 

z — w 

+ ... 


for ji / j 2 and si, S 2 = ±1. Moreover, when a • {Eej) = 
1, i.e. Q: zb Gjf G AgQ('2r-|-l) 5 

E^(z)E^^^ (w) = —b£;«±e, (^U)) + ... 


where e(m, n) = (—1)"^*^ is defined in (BIO). 


The (normal ordered) energy-mo mentu m tensor can 
be turned from the Sugawara form (2.12) to the usual 
bosonic and fermionic one 


T{z) 


1 


J2h\z)H\z) + J 2 E°{z)E-'^{z) 

_i=l ckGA 


2{N - 1) 

r 

+ E {z)E-^i (z) + E-^i {z)E^i (z) 


j=l 


= -Y(l>{z) ■ d4>{z) - Y^^+Hz)d^^^+Hz). (C4) 

There are only two non-trivial primary fields Y nnd a. 
The fermion sector Y consists of the 2r Dirac fermions 


in (B23) as well as the remaining Majorana 
fermion '/^’^+^TThe a sector consists of 2'^ twist fields 




(C5) 


= span < exp i 


i E ^ 4>^ I : Sj = 0,1 


i=i 


which represents so{2r + l)i according to the spinor rep¬ 
resentation. Their conformal dimensions are given by 
= 1/2 and ha- = (2r -b 1)/16. 


Appendix D: Ze parafermion model 

Here we represent the Ze parafermions using bosonized 
fields and Majorana fermions in the so(9)i CFT. We fo¬ 
cus on a single Majorana wire containing 9 right moving 
real fermions. The CFT is fractionalized using the con¬ 
formal embedd ing into so(9)i D so(3)/ x so(3)/ (see 
section 


III Bil l. Each 50(3)3 sector is then further de¬ 
composed into 50(2)3 X “Ze” using th e coset construc¬ 


tion “Ze” = 50 ( 3 ) 3 / 50 ( 2)3 (see section IIIB2). We now 
provide a more detail description of the Ze parafermion 
sector. We will focus on the one in 5o(3)/. 

First we pair six Majorana channels into three Dirac 
fermions and bosonize = (7/^ -b i '/^)/\/2 = , Y = 

(7/^ -b iY^)/\/2 = and -b iY^)/\/2 = 

The Lagrangian density of the boson fields are given in 
(3.68). Like the 50(A) 1 case, extra care is required so 
that the Dirac fermions satisfies the appropriate mu¬ 
tual anticommutation relations. Here we use a slightly 
different but more convenient convention 


0/z)0'^ (rc)y = —6'^^ fog(^ — w) 

0 if ^ = j 
1 if i — j = 1 mod 3 
— 1 if i — j = —1 mod 3 


(Dl) 


5*^' = 


SO that the constant phases have a threefold cyclic 
symmetry. The 50 ( 2)3 sub-theory is generated by the 
“charged” boson (j)p = {cjY -\r (jY 0^)/3. It satisfies 


{^p{z)(t>p{w)) = -- log( 2 : - w). 


(D2) 


The remaining “neutral” bosons —4>p are linearly 

dependent = 0 and obey the OPE 

{<i>a{zWa(w)) = - - Y log(^ -W) + 

(D3) 

The “charge” and “neutral” sector completely decoupled 
so that {(j)p{z)(j)i^{w)) = 0. Lastly, there are three remain¬ 
ing Majoranan fermions in the 5o(9)i theory. They 
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completely decouple with 4>a- and 0^. Although the ver¬ 
tex anticommutes with this has no effect on 

any of our derivations. More importantly the “neutral” 

vertices commute with the remaining fermions^ _ 

In section IIIB 2, we defined the Ze parafermion ( 3.75| ) 


^ (D4) 


which is part of the so{3)^ current (see ( 3.74[ )). It gen¬ 
erates the rest of the Zq parafermions 


= 


^/lE 

‘li 


i=i 


- cos {(1)1 - 4)1) 4)^ (D5) 

-COs(</.^-</.3)y-COs(<^3_^l)^8] 

^ g/® = (g'i)^, g/° = g/® = i 

where and ^p^^^ = ^p^^p^^p^. Their conformal 

dimensions 

m (6 — m) 


fl^m — 

as well as the fusion rules 

^r, 

^^{z)^^'iw) = 


{z - 

1 

^ 2/l^T77. 

2h^rj 


(D 6 ) 


^m+m'(y,)+ __ (D7) 


■{z - w)^Tze + . . . 


1 


match with the known result by Zamolodchikov and 
Fatee^^^, for T^g the energy-momentum tensor ( |3.73| ) 
with central charge czq =5/4 and 


' (m + m ')!(6 — m )!(6 — m')\ 
m\m'\{6 — m — m')! 6 ! 


(D 8 ) 


which are identical to the modular 5'-matri:jP^ of the 
Qn WZW CFT. The fusion matrices that charac¬ 
terize fusion rules a x b = turned be 

determined by ^'-matrix throught the Verlinde formulsP^ 

(|3^ 


lyrs _ •^Sis'Ss2S'Sss' 
^^SiS2 ~ 2-^ , 

, *^0s' 


(E3) 


The Gn state is Abelian and carries four anyon 
types l,' 0 ,s+, 5 _ when A" is a multiple of four. It is 
non-Abelian otherwise and carries three anyon types 
l,'0,cr when A is 2 mod 4, or seven anyon types 
1, (a+ 7 +,/3, 7 _, (a_, / when A is odd. The quasiparticle 
exchange statistics and quantum dimensions are 
summarized in t able [TT| and [ml The total quantum di¬ 
mensions V = \/X]x given by 


= 


2 for A even 

2 csc( 7 r/ 8 ) for A odd 


(E4) 


where csc(7r/8) = ^/4 + 2v^. 

The ^'-matrices of Gat for A = 2r even ar e well- known 
and are given by those of the SO{r)i states 


^Gn = 
^Gn = 


1 


'^Gn 

1 

'^Gn 


111 1 
11 - 1-1 
1 -1 
1 -1 

1 1 V 2 
1 1 -V2 
V 2 -\/2 0 


for A = 4n, (E5) 

for A = 4n + 2 . (E 6 ) 


The ^'-matrices for the odd A ca ses are modification of 
the Gg = SO(3)s prototype ( |3.59[ ) 




7r(2si + 1)(252 + 1) 


(E7) 


where Sj = 0,1/2,1, 3/2, 2, 5/2, 3 label the seven anyon 
types 1, (a+, 7 +, /3, 7 _, a_, / (see table Eor A = 9 + 2r 
mod 32, the ^'-matrix of Gat is given by 


Appendix E: The AS-matrices of the Gn state 

The surface topological orders of the time reversal sym¬ 
metric gapped coupled wire model are described in sec¬ 
tion El Ther e ar e thir ty tw o distinct topological states 
defined in eq.(|4.6|) and (|4.9|), which we repeat here. 


for A = 2r 


G = / ('El') 

^ \ S' 0 ( 3 ) 3 KlbS' 0 (r)i, for Af = 9 + 2r ■ ^ ^ 

In this appendix we summarize the modular properties 
of these sta tes. In particular we present there braiding 
A^-matrices (4.1) 

= (E2) 


SG^=r'S%\r/2])T- 


(E 8 ) 


where \r/2] > r j2 is the smallest integral ceiling of r/ 2 , 
S^{n) is the ^S-matrix when r = 2 n is even 




(E9) 


and T is the operator that flips the fermion parity of 
a- and 7 + 7 _ 


J^ = 




'v ^ 


1 \ 


(ElO) 
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